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ABSTRACT 

We present an unbiased census of deeply embedded protostars in Perseus, Serpens, and Ophi- 
uchus, assembled by combining large-scale 1.1 mm Bolocam continuum and Spitzer Legacy surveys. 
We identify protostellar candidat es based on their m i d-infrared propertie s, cor relate their p o sition s 
with 1.1 mm core positions from lEnoch et all ()2006f ). lYoung et all (|2006f ). and lEnoch et al.l (|2007| ). 
and construct well-sampled SEDs using our extensive wavelength coverage (A = 1.25 — 1100 ^m). 
Source classification based on the bolometric temperature yields a total of 39 Class and 89 Class I 
sources in the three cloud sample. We compare to protostellar evolutionary models using the bolo- 
metric temperature-luminosity diagram, finding a population of low luminosity Class I sources that 
are inconsistent with constant or monotonically decreasing mass accretion rates. This result argues 
strongly for episodic accretion during the Class I phase, with more than 50% of sources in a "sub-Shu" 
[dM/dt < 10^^ Mq yr~^) accretion state. Average spectra are compared to protostellar radiative 
transfer models, which match the observed spectra fairly well in Stage 0, but predict too much near-IR 
and too little mid-IR flux in Stage I. Finally, the relative number of Class and Class I sources are 
used to estimate the lifetime of the Class phase; the three cloud average yields a Class lifetime 
of 1.7 ± 0.3 X 10^ yr, ruling out an extremely rapid early accretion phase. Correcting photometry 
for extinction results in a somewhat shorter lifetime (1.1 x 10^ yr). In Ophiuchus, however, we find 
very few Class sources (-/Vciass oZ-^ciass i ~ 0.1 — 0.2), similar to previous studies of that cloud. The 
observations suggest a consistent picture of nearly constant average accretion rate through the entire 
embedded phase, with accretion becoming episodic by at least the Class I stage, and possibly earlier. 
Subject headings: stars: formation — ISM: clouds — ISM: individual (Perseus, Serpens, Ophiuchus) 
- submillimeter - infrared: ISM 
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1. INTRODUCTION 

The problem of how low mass stars like the sun form 
has been studied extensively over the last few decades. 
Compared to more evolved protostars and pre-main se- 
quence objects, however, the earliest stages of the star 
formation process, from the formation of dense cores 
through the main mass accretion phase, are relatively 
poorly understood. This lack of information is due in 
large part to the difficulty of observing young, deeply 
embedded sources, which are shrouded within dense pro- 
tostellar envelopes and only observable via reprocessed 
emission at mid-infrared to millimeter wavelengths. Fur- 
thermore, most previous observations of deeply embed- 
ded objects have naturally focused on a small number 
of very bright or well known sources, due to the sensi- 
tivity and resolution limitations of long-wavelength sur- 
veys. Understanding the formation of typical stars re- 
quires complete samples of young objects, over molecular 
cloud scales. 

Currently, the details of the early evolution of proto- 
stellar sources are extremely uncertain, including mass 
accretion rates during the Class and Class I phases. 
In addition, measurements of the timescales associated 
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with the earliest stages vary con siderably, ranging from 
lOfto 10^ yr for prestellar cores (|Ward-Thompson et al.l 
I2007D and from lOMo a few xlO^ yr fo r Class 
(| Andre fc Montmerlel[T99l IVisser et al.ll2002D . In fact, 
the association of Class and Class I with distinct 
evolutionary stages is sti ll a matter of debate (e.g., 
iJavawardhana et all l2001f ). Large surveys at mid- 
infrared to millimeter wavelengths, where the SEDs of 
embedded sources peak, are essential for understanding 
how protostars evolve through their earliest stages. In 
addition to providing complete samples of young sources, 
large surveys are also important for characterizing vari- 
ations in the star formation process wit h environment. 

Wi t h a few notable recent exceptions (|j0rgensen et al.l 
l2007t iHatchell et al.l l2007al) . previous samples of very 
young protostars have typically been compiled from 
many different surveys, and suffered from systematics, 
unquantified environmental effects, and small number 
statistics. We recently completed large-scale 1.1 mm 
continuum surveys of Perseus, Serpens, and Ophiuchus 
with Bolocam at the Caltech Submillimeter Observatory 
(CSO). Maps have a resolution of 31" and cover 7.5 deg^ 
in Perseus (140 pc^ at our adopted cloud distance of 
d = 250 pc), 10.8 deg2 in Ophiuchus (50 pc^ at d = 
125 pc), a nd 1.5 de g^ in Serp ens (30 pc^ at d = 260 pc 
(lEnoch et al. 2006; Young et al. 2006: lEnoch et al.ll200 
hereafter Papers I, II, and III, respectively). These Bolo- 
cam surveys complement large Spitzer Space Telescope 
IRAC and MIPS maps of the same clouds from the 
"From Molecular Cores to Planet-forming D isks" Spitzer 
Lega cy program ("Cores to Disks" or c2d; lEvans et all 
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Millimeter emission traces the dust in dense starless 
cores and protostellar envelopes, and provides a measure 
of core and envelope properties, including sizes, masses, 
and spatial distribution. Spitzer IRAC and MIPS ob- 
servations are complementary in that they provide in- 
formation about the properties of any young protostars 
embedded within dense cores. Combining these data en- 
ables us to assemble a mass limited, unbiased census of 
the youngest star-forming objects in three different envi- 
ronments, including prestellar cores. Class 0, and Class I 
protostars. 

In Paper III we looked at how the global cloud envi- 
ronment in fluences the properties of star-forming cores 
( Enoch et al. 2007). In a companion paper to this work 
(|Enoch et al. 2008), we examine the properties of prestel- 
lar and protostellar cores in Perseus, Serpens, and Ophi- 
uchus, focusing on the prestellar core mass distribution 
and the lifetime of prestellar cores. A similar analysis 
comparing the c2d Spitzer data to large SCUBA maps 
of Perseus and Oph i uchus has recently been carried out 
bv |j0rgensen et all (|2007l [2008), focusing on the differ- 
ence between starless and protostellar cores, as well as 
cloud properties such as the star formation efficiency 
and how it varies with s patial clustering. We follow 
iDi France sco et all (|200 7^ in defining millimeter cores 
containing a compact luminous internal source (i.e., an 
embedded protostar) as "protostellar cores", regardless 
of whether the final object will be stellar or sub-stellar 
in nature. We use "starless cores" to refer to dense cores 
without an internal luminosity source, and "prestellar 
cores" as starle ss cores that are li kely to be gravitation- 
ally bound fsee lEnoch et al.ll2008f) . 

In this work we exploit the combined power of millime- 
ter and mid- to far-infrared observations to study the 
evolution of young protostars embedded within the pro- 
tostellar cores. Extensive wavelength coverage from A = 
1.25 — 1100 fim allows us to trace the evolution of proto- 
stellar sources in their main mass accretion phase, from 
format ion through the e nd of the e mbedded pha s e. We 
follow iRobitaille eFaLl (|2006f ) and ICrapsi et all (|2008l ) 
in using "Stage" (e.g.. Stage 0, I, II, III) to refer to 
a source's true physical nature, regardless of its ob- 
served properties, while the corresponding "Class" refers 
to the observational classification, typically based on 
the near- to mid-in frared slope aiR (jAdams et al.lll987l : 
iGreene et all 119941) or on the bolometric temperature 
(|Mvers fc Laddlll993h . 

We assume that Stage and Stage I refer to an evo- 
lutionary sequence of embedded protosta rs with Mh, < 
Menv and > Me„t,, respectively (e.g. lAndrlll994f ). 
where is the protostar mass and M^nv the enve- 
lope mass. Similarly, we assume Stage II refers to pre- 
main sequence st ars with very little remaining envelope 
{M^nv < O.IM0; ICrapsi et al1l2008f l. Our adopted def- 
initions of various Classes and Stages are summarized 
in Table [TJ Ideally, Class I wo uld directly corr e spond 
to S tage I, etc., b ut as shown by Robitaillc et al. (2006) 
and ICrapsi et al.l (pOOSi ). geometric effects can cause, for 
example, Stage II sources to be classified as Class I. 

In ^J2]we describe the 1.1 mm and Spitzer infrared (IR) 
data, including how protostellar candidates are identi- 
fied and their association with 1.1 mm cores. We cal- 
culate bolometric luminosities and temperatures, enve- 
lope masses, and discuss completeness (^. Protostel- 



lar classification methods are compared in Sj4| Spectral 
characteristics of Class 0, I and II sources detected at 
1.1 mm are discussed in ^|S1 including selected individual 
sources f ij5.2p . In Sj6]we explore alternative classification 
schemes made possible by our large sample, and com- 
pare the average observed spectra to protostellar models 
in SjTj Sources are placed on a bolometric temperature- 
luminosity diagram for comparison to protostellar evolu- 
tionary models and to study their luminosity evolution, 
mass accretion rates, and envelope evolution (SjSj). Fi- 
nally, in ^j9]we calculate lifetimes for the Class phase. 

2. COMBINING BOLOCAM AND SPITZER C2D DATA 

Both Bolocam 1.1 mm and Spitzer maps were designed 
to cover down to a visual extinction of Ay ^ 2 mag in 
Perseus, Ay > 3 mag in Op hiuchus, and Ay ^ 6 mag in 
Serpens (|Evans et al.l 120031 ). The actual overlap in area 
between Bolocam and IRAC maps is shown in Figure 1 
of Papers I, II, and III for Perseus, Ophiuchus, and Ser- 
pens, respectively. Catalogs listing c2d measured Spitzer 
fluxes of all sources in each of the three clouds, as well as 
near-infrared fluxes for sources in the 2MA SS catalogs, 
are a vailable through the Spitzer database (|Evans et al.l 
|2007( ). We utilize wavelength coverage from A = 1.25 to 
1100 fim, using 2MASS (1.25, 1.65, 2.17 ^im), IRAC (3.6, 
4.5, 5.8, 8.0 fim), MIPS (24, 70, 160 ^m), and Bolocam 
(1.1 mm) data. Note that 160 fim flux measurements are 
not included in the c2d delivery catalogs due to substan- 
tial uncertainties and incompleteness, but are included 
here when possible. Photo metry at 160um is d iscussed 
in RebuU et al. (200'3) and lHarvev et al.l (j2007bl ). where 
160/im fluxes for point sources in Perseus and Serpens 
are also given. 

Basic data papers describe the processing and analysis 
of the Spitzer IRAC and MIPS maps of Perseus, Ser- 
pens, and Ophiuchus, as well as presenting general prop- 
erties of the sources in each clo ud, such as colo r -color 
and color-ma gnitud e diagrams (iJorgensen et al.] 120061: 
Harvey et al! 120061 : iRebull et all 120071 : iHarvev et all 
2007bf) . The young stellar objec t (YSO) popu lation in 
Serpens is discussed in detail bv IHarvev et al.l ([2P07a). 
Here we are most interested in the young protostellar 
sources that are most likely to be embedded in the mil- 
limeter cores detected with Bolocam. For the following 
we will use the term "embedded protostar candidate" in 
general to encompass candidate Stage I and younger ob- 
jects. Although it will become apparent in i i2.1l that our 
criteria also pick up a number of Stage II sources, in gen- 
eral we focus on sources with evidence for a protostellar 
envelope (Stage and I). 

Figure [T] shows the result of combining Spitzer and 
Bolocam data, for a few examples of embedded protostel- 
lar sources in each cloud. Images are three-color Spitzer 
maps (8.0, 24, 70 /im), with 1.1 mm contours overlaid. 
The Bolocam ID of the associated 1.1 mm core is given 
in each panel, as well as the embedded source IDs from 
Tablesllll We use the 2MASS, IRAC, AilPS, and Bolo- 
cam data to construct complete SEDs from 2 to 1100 /im 
for each candidate protostellar source, shown in the lower 
panels (see also . Additional SHARC II 35 nm fluxes 
jmet al.1 12001 and SCUBA 850 fim fluxes (jKirk et all 
I2OO6I open circles) are included when available. Modified 
black-body {Si, cx vl^By{T)) curves for a temperature of 
15 K and (3=1 are shown for reference in Figure [T] (dot- 
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Fig. 1. — Three-color Spitzer images 8 fim (blue), 24 (green), and 70 fim (red) of selected embedded protostar candidates in Perseus, 
Serpens, and Ophiuchus, with Bolocam 1.1 mm contours. The Bolocam ID of the centered core and embedded protostar IDs from Tables[2}j4] 
given, and contour intervals for the 1.1 mm emission are 3, 6. ..15, 20, 25, ...o". Spectral energy distributions (SEDs) of the candidates are 
plotted in the lower panels. SEDs include 2MASS, IRAC, and MIPS photometry from the c2d database and 1.1 mm fluxes (asterisks). 
Published 350 fim and 850 /im points (open circles) are also shown when available. A modified blackbody spectrum (T = 15 K, /3 = 1) is 
overlaid for reference. 
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ted lines). 

Protostellar sources in our sample may be isolated 
(e.g., Per-Bolo 57) or lie in crowded regions (e.g., Per- 
Bolo 49). Approximately 20%-50% of the time, more 
than one protostellar source lies within a single 1.1 mm 
core (20% of the embedded protostar sample in Ophi- 
uchus, 40% in Perseus, and 55% in Serpens). Much of 
this difference is certainly due to the lower resolution of 
Bolocam (31") compared to Spitzer (7" at 24 /im), caus- 
ing nearby Spitzer sources to be blended in the 1.1 mm 
map, although the envelopes detected at 1.1 mm are also 
physically more extended than the region emitting at 
Spitzer wavelengths. Sometimes the SEDs of such mul- 
tiple sources look similar to each other (Per-emb 18 and 
22 in the lower-left panel of Figure [1]) , and sometimes 
quite different (Per-emb 49 of the same panel). 

Although the coverage of the IRAC, MIPS and Bolo- 
cam maps overlaps nearly perfectly for our purposes, 
there are a few cases in which embedded protostar candi- 
dates are outside the boundaries of the MIPS 70 /im map 
(three sources in Serpens) or 1.1 mm map (one source in 
Perseus) . In addition, a number of bright sources in each 
cloud are saturated in the 24/^m or 70/j,m c2d maps. In 
these cases we substitute IRAS 25/im or 60/im fluxes is 
the source does not appear to be blended in a visual in- 
spection of the the IRAS maps. 

The 160 /xm maps are often saturated near bright 
sources and in regions of bright extended emission, such 
as near bright clusters of sources. Reliable 160 /zm fluxes 
are especially difficult to determine in crowded regions, 
due both to the large beam size (40") and to saturation 
issues. The lack of 160 /im data is most problematic in 
Ophiuchus, where the 160 maps are saturated in all 
of the dense source regions. Even for isolated sources, 
the measured 160 /im flux density, determined from a 
point spread function (PSF) fit, may be underestimated 
if the source is extended. The effects of these issues on 
our analysis are discussed in more detail in the appendix 

(§©• 

2.1. Identifying Embedded Protostar Candidates 

We form a sample of candidate embedded protostars 
from the c2d catalogs; the first cut is based on the source 
"class." All sources in the c2d catalogs are assigned a 
class parameter based on colors, magnitude s, and stellar 
SEP fits ( see the c2d Delivery Document () Evans et al.l 
|2007() and iHarvev et al.l ()2007aD '). Class parameters in- 
clude "star", "star-hdisk" , "YSOc" (young stellar object 
candidate), "red", "rising", "Gale" (galaxy candidate), 
etc. Embedded protostars will generally be a subset of 
YSOc sources, but some of the most embedded may also 
be assigned to the "red" class if they are not detected in 
aU IRAC bands. 

Thus, we begin by selecting all sources from the c2d 
database that are classified as "YSOc" or "red." From 
this list, we keep sources that meet all of the following 
criteria: 

(a) flux density at 24 /im (524/im) > 3 mJy 

(b) S'24/im > 5a/fl + 8 mJy, where am is the near- 
to mid-infrared spectral index, determined by a least- 
squares fit to photometry between 2 and 24 /xm. This 
criteria is motivated by a co mparison to the care fully 
vetted Serpens YSO sample of iHarvev et al.l (|2007aD . 

(c) i^S'24/im > i^Ss^ui, i-e., the SED is rising from 8 to 



24 fj,m in vSi, space. 

(d) S'24/im must be of high quality, i.e., signal to noise 
(S/N) greater than 7. 

(e) S'24/im is not a "band-filled" flux. For sources not 
originally detected in all Spitzer bands, a flux or up- 
per limit is measured at the source position (band-flUing, 
IHarvev et al.ll2007al ). Because the resolution is lower at 
24 /im than at the shorter wavelengths, some IRAC-only 
sources have unreliable band-flUed fluxes at 24 /im (e.g., 
sources are confused with the PSF wings of a nearby 
source). 

In addition to sources that meet the above criteria, we 
include any 70 /im point sources not classified as galaxy 
candidates ("Gale"). Note that these 70/im sources need 
not be classified as "red" or "YSOc" . In each cloud, a 
number of deeply embedded sources that are bright at 
70 /im but very weak at 24 /xm (e.g., HH211 in Perseus) 
are recovered by this last criteria (5 in Perseus, 4 in 
Serpens, and 3 in Ophiuchus), as are a few very bright 
sources that are saturated at 24 /im (6 in Perseus, 2 in 
Serpens, and 6 in Ophiuchus; these are often classified as 
"rising"). 
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Fig. 2. — Plot of S'24/im versus spectral index am for sources in 
Serpens, to demonstrate selection criteria for the embedded pro- 
tostar candidate samples. Plus symbols indicate the sam- 
ple from which the candidate embedded protostars is drawn, in- 
cluding all catalog sources labeled as "YSOc" or "red" as well as 
non-"Galc" 70/xm sources. Lines show 24 fim flux cuts imposed 
by criteria (a) and (b), which remove the majority of spurious or 
background sources. Boxes indicate the final embedded protostar 
sample, after applying criteria (a)-(d). Smal ler blue di amonds show 
the carefully vetted YSO sample from Ha rvey et al.l jgOOS). The 
two samples agree fairly well, but because we select for embed- 
ded sources we miss many of the more evolved YSOs in the vetted 
YSO list (Class II/III). There are a few sources in the embedded 
protostar sample not in the Ifarvey et al. (2006) list, two of which 
are likely to be real (see text), while the others are rejected when 
examined by eye. 
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Figure [2] plots S'24/im versus am for embedded candi- 
dates in Serpens, where "+" symbols indicate the orig- 
inal ("YSOc"-|-"red"-|- 70/im sources) sample for Ser- 
pens, and boxes indicate our embedded protostar can- 
didates after applying criteria (a)-(e). For compar- 
ison, diamond s indicate the carefu lly vetted Serpens 
YSO list from iHarvev et all (|2007a[ ). The majority of 
"-|-"-only poi n ts, wh ich were rejected as true YSOs by 
IHarvev et al.l ()2007aD . are removed by criteria (a) and 
(b) (shown as solid lines). More than half of the vetted 
YSOs do not appear in the candidate embedded proto- 
star sample (diamond, no box); the majority of these are 
rejected by the rising SED criteria (c) and as they are 
primarily classified as "star-|-disk" we do not expected 
them to be embedded. There are a few sources in our 
sample that are not in the vetted YSO sample (box, no 
diamond). A few of these were easily identified as non- 
protostellar, and rejected, when examining by eye. We 
do identify three embedded protostar candidates in Ser- 
pens that are not in the vetted YSO list but seem to be 
associated with 1.1 mm emission (see WZ.'Z^ : one "red", 
and two "rising" . 

We include red sources to ensure that we identify the 
most embedded protostars, but we need to reject the 
large percentage of these sources that are likely to be 
galaxies. As shown in Figure[ni criteria (a) and (b) are ef- 
ficient at eliminating extragalactic sources. Here S'24^m 
is plotted versus am for all "Gale" galaxy candidates in 
Serpens. As in Figure [H lines show criteria (a) and (b) 
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Fig. 3. — Plot of S24iim versus spectral index aju for galaxy 
candidates ("Gale") in Serpens. Selection criteria (a) and (b) are 
overlaid, as in Figure [2] These selection criteria eliminate most of 
the parameter space inhabited by extragalactic sources. 

After forming our embedded candidate samples based 
on the above criteria, the images and SEDs of each source 
are examined by eye to remove any galaxies that are 
extended in the near-infrared, and other obviously non- 
embedded sources. SEDs include 2MASS, IRAC, MIPS 
and 1.1 mm fluxes. In some cases, there is no available 
point-source flux at 70 or 160 /im even if there is emis- 



sion at the position of the source, generally because the 
source is extended at these wavelengths. In these cases a 
flux density is measured by hand, if possible, using aper- 
ture photometry. Large uncertainties (50% or more) are 
associated with these band-filled flux measurements. 

2.2. Association with a 1.1 mm Core 

The next step after assembling a sample of candidate 
embedded protostars is to determine which are associ- 
ated with 1.1 mm emission. The correlation between 
ca ndidate protostar a nd 1.1 mm core positions was done 
in lEnoch etlHI (l2008l see ?2.3 and Figure 2), following 
a similar analvsis bv lj0rgensen "eFall (I2OOI : we summa- 
rize the results here. We found that the coldest proto- 
stellar candidates (Tt,oi < 300 K, see are primarily 

located within 1.0 x Oimm of a millimeter core position, 
where Oimm is the angular FWHM size of the 1.1 mm 
core, as expected if they are deeply embedded. Based on 
that analysis, a given embedded protostar candidate is 
assumed to be associated with a millimeter core if it is 
located within 1.0 x 9imm of the 1.1 mm core centroid 
position. 

We also found in lEnoch et al.l (|2008f) , based on anal- 
ysis of a spatially random distribution of sources, that 
we can expect approximately 5 false associations with 
1.1 mm cores in each cloud. Using a more restrictive cri- 
teria (e.g., 0.5 X Oimm) would reduce the number of false 
associations, but would likely miss at least a few embed- 
ded protostars. For reference, adopting 0.5 x Oimm would 
result in 10 fewer embedded protostars in Perseus, 7 in 
Serpens, and 16 in Ophiuchus. 

If an embedded protostar candidate is found to be as- 
sociated with a millimeter core by the above criteria, the 
1.1 mm flux is included in the protostellar SED, and is 
used to calculate an envelope mass {Menv, ^3.2p . If a 
candidate is not located within 1.0 x Oimm of a millime- 
ter core position, we calculate a flux or upper limit from 
the original 1.1 mm map (although sources with upper 
limits at 1.1 mm are not included in the flnal lists of 
embedded protostars given in Tables [2]-[4]) . 

In a number of cases (10 in Perseus, 3 in Serpens, and 
11 in Ophiuchus) there is clearly 1.1 mm flux at the pro- 
tostar position that was not identified as a core in the 
original 1.1 mm source extraction because it is below ha 
(where cr is the local rms noise, see Paper I) or is in a con- 
fused region of the map. If this emission exceeds 3cr, we 
measure a "band-filled" 1.1 mm fiux using a small aper- 
ture (30"-40"). We also re-compute the 1.1 mm flux 
density in small apertures for all sources in regions of 
blended 1.1 mm emission. When more than one proto- 
star candidate is associated with a single compact mil- 
limeter core, we divide the total 1.1 mm flux of the core 
equally between the protostellar sources. 

The approach described above means that the 1.1 mm 
flux is not measured in a fixed aperture for every source. 
In crowded regions or for sources with close neighbors, 
the 1.1 mm flux is measured in a 30"-40"diameter aper- 
ture centered on the Spitzer source position. These 
sources are indicated by a footnote in Tables [2H11 For 
other sources, the Bolocam 1.1 mm core flux from Papers 
I-III is used. Typically, but not always, the total core 
flux is used here (integrated in the largest aperture, from 
30" — 120" diameters in steps of 10", that is smaller than 
the distance to the nearest neighboring source); fluxes 
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in 30, 40, 80, and 120" apertures can also be found in 
Papers I-III. 

Approximately 50% of the embedded protostar candi- 
dates in each cloud are lacking 1.1 mm emission, even 
after re-examining the 1.1 mm maps at each source po- 
sition. Most of these sources appear to be Class II ob- 
jects, with little or no remaining envelope. Any 1.1 mm 
emission from these objects, therefore, is likely below our 
detection limit of ~ 0.1 Mq (§ 13. 3|) . A few of these may 
actually be low luminosity sources that are truly embed- 
ded, with low mass envelopes below our detection limit. 
For this reason, we only claim completeness to embed- 
ded protostars with M^,n„ > 0.1 M f^. Note that this 

corresponds well to the ICrapsi et al.1 (|2008D definition of 
Stage I. 

3. PROPERTIES OF EMBEDDED PROTOSTARS 

3.1. Bolometric Luminosity and Temperature 

Our extensive wavelength coverage allows us to con- 
struct well-sampled SEDs for the embedded protostar 
candidates in all three clouds, from which we calculate 
a bolometric luminosity (Ltoi) and bolometric tempera- 
ture (Thoi) for each source. The bolometric luminosity is 
calculated by integrating the SED {S^) over frequency: 



L 



bol 



(1) 



The bolometric temperature is defined as the tempera- 
ture of a blackbody with the same mean frequency as the 
source SED, and is given by 

Tboi = 1.25 X 10-"(i^) K (2) 

(jMvers &: Laddl[l993[ ). where the mean frequency is 



(3) 



J S^dv 

Two methods for approximating the integrations over 
frequency for finitely sampled SEDs (midpoint and pris- 
moidal) are discussed in the appendix 

Tables [2]-[4] list the bolometric temperatures and lu- 
minosities derived for all embedded protostars in each 
cloud. As we are primarily interested in young, em- 
bedded objects, only sources with detectable 1.1 mm 
emission, which likely have substantial envelopes, are 
included. Sources are listed by increasing T^oi, and 
identified as, e.g., "Per-emb#", as well as by their c2d 
name (SSTc2dJ...), which also gives the position. When 
sources are saturated at 24 or TO/im, the IRAS 25 or 
60/Ltm flux is utilized when not affected by blending. In 
these cases a note is made in Tables [2H4l The correction 
of saturated fluxes can increase the luminosity by more 
than a factor of two; for saturated sources with no reli- 
able IRAS flux (also noted in TablesUHl]), the luminosity 
will be an underestimate. 

All Tiioi and L^oi values quoted use the midpoint in- 
tegration method. The difference between the values 
calculated by the midpoint and prismoidal integration 
methods (given in parentheses in Tables [2H4|) gives a 
more realistic measure of the uncertainties in T\,oi and 
Liioi than the formal fitting errors (which are typically 
10%). Uncertainties can be larger than 50% depending 
on whether or not a 160 /im flux is available, and there 
is an additional systematic uncertainty of 15-25% from 
finite sampling errors 



Systematic errors introduced by missing 160 /^m fluxes 
are investigated in ^A) Without a 160 /xm measurement, 
Tiioi will almost certainly be an overestimate for very cold 
sources, which may bias our classification of protostellar 
candidates (Sj4]). Ophiuchus will be most affected, as only 
four sources in that cloud have reliable 160 ^va fluxes. 

3.1.1. Correcting for Extinction 

One might argue that we should correct our photom- 
etry for extinction before calculating Li^oi and Tboz- The 
effects of extinction are typically ignored for Class and 
Class I sources, for which one might expect the fore- 
ground extinction to be negligible relative to the ef- 
fect of the envelope itself. While these observed val- 
ues are m ost easily compare d to the majority of previ- 
ous work. IChen et al] (|1995( l. who connected Tboj to the 
classes defined by a, did correct the observed flux den- 
sities for extinction before computing T^o;. Therefore, 
we compute exti nction corrected values for comparison 
(see lEvans et al.il2008l for a more complete discussion of 
the effects of dereddening) . We use the mean extinction 
to all Class II sources in each cloud {Ay = 5.92 mag in 
Perseus, 9. 57 mag in Serpens, an 9.76 mag in Ophiuchus; 
lEvans etl il. 2008) to deredden the photometry of Class 
and Class I sources. 

As dereddening has a greater effect on shorter wave- 
length photometry, correcting for extinction tends to in- 
crease both Tfjoj and Lboi- Throughout this paper we 
will primarily use observed values, and not extinction 
corrected values, but the effects on derived values such 
as the Class lifetime (§ [9]) will be noted. When used, 
the extinction corrected bolometric luminosity and tem- 
perature are indicated by T^'^; and L'^^^, respectively. 



3.2. Envelope Mass 

The envelope mass of candidate protostars, Mg. 
calculated from the flux density at 1.1 mm, Simm' 



M 



d^Sir 



IS 



(4) 



where d is the cloud distance, Kimm — 0.0114 cm^ g~^ 
is the dust opacity per gram of gas at 1.1 mm, and 
Birnm IS the Planck function at a dust temperature of 
T/j. The opacity is interpolated from Table 1 column 5 
of Ossenkopf & Hcnning ( 1994) for dust grains with thin 
ice mantles. 

We assume a dust temperature of To = 15 K for proto- 
stellar envelopes, consistent with average isothermal dust 
temperatures found from radiative transfer mod els of a 
sample of Class and_Class I protostars (Shirley " et al.l 
120021: lYoung et al.ll2003l) . The isothermal dust tempera- 
ture is the temperature that, when used in an isothermal 
mass equation (e.g., Eq. (|4]), above) yields the same mass 
as does a detailed radiative transfer model that accounts 
for temperature gradients. Dust temperatures will be 
higher than 15 K close to the protostar, but the major- 
ity of the envelope mass is in the outer, cooler, regions of 
the envelope. A dust temperature of 10 K would result in 
an increase in masses by a factor of 1.9, while To = 20 K 
would decrease masses by a factor of 1.5. 

Envelope masses of embedded protostars are listed in 
Tables [5H11 If the source is associated with a distinct 
core, the Bolocam identification from Papers I-III is 



Protostars in Serpens, Perseus, and Ophiuchus 



7 




Fig. 4. — Near- to mid-IR spectral index (aju) versus bolometric temperature (Tjo;) for the candidate embedded protostar samples in 
Perseus, Serpens, and Ophiuchus. Bold symbols indicate sources that are associated with 1.1 mm emission, and thin symbols denote those 
with upper limits at 1.1 mm. Standard class divisions for both Ttoj and am are shown. The two methods agree fairly well for Class II and 
"warmer" Class I sources, but very cold (Ti,o; < 100 K) sources have a large range of am values. 



given in the last column. If the source is not associated 
with one of the originally identified cores, but rather the 
flux has been "band-filled" at 1.1 mm, then no Bolocam 
ID is given. In these "band-filled" cases, the measured 
1.1 mm flux can be easily re-constructed from the mass. 

3.3. Completeness 

Because we require a detection at 1.1 mm to be in- 
cluded in the final source lists, we are clearly incom- 
plete to objects more evolved than Stage I, which do 
not have a substantial protostellar envelope. We do de- 
tect a few Class II sources at 1.1 mm; these are dis- 
cussed in i j5.ll The point-source detection limits of our 
1.1 mm surveys limit our sensitivity to Class I sources 
with Menv > 0.09 Mq in Perseus, Me„„ > 0.07 Mq in 
Serpens, and Alenv ^ 0.04 Mq in Ophiuchus. These are 
5cr point source detection limits; for very extended en- 
velopes the completeness limits will be higher. For sim- 
plicity we take a detection limit of O.IMq for all clouds. 

As the 1.1 mm detection requirement is more restric- 
tive than our 24 /itm flux criteria, we explore the possibil- 
ity that we are missing some low luminosity embedded 
sources that are below our 1.1 mm detection threshold. 

Taking the 1.1 mm 5ct detection limits for each cloud 
(75 mJy in Perseus, 50 mJy in Serpens, and 110 mJy in 
Ophiuchus), and assuming the spectrum of a modified 
blackbody, = vl^B^iTo), with Tq = 15 K, we can 
estimate the minimum detectable bolometric luminosity 
in each cloud. For P — 1, the minimum Ltoi is 0.02 Lq in 
Perseus, 0.01 L0 in Serpens, and 0.01 Lq in Ophiuchus. 



Assuming Tjj = 20 K lowers these values by approxi- 
mately a factor of 2, while taking (3 — 2 increases them 
by a factor of 4. Although these are very rough estimates, 
they agree fairly well with the lowest observed bolometric 
luminosities for sources with 1.1 mm detections (0.04 Lq 
in Pers eus, 0.05 Lq in Serpe ns, and 0.01 Lq in Ophi - 
uchus). iDunham et aH (|2008[ ) and lHarvev et al.l (|2007bl ) 
demonstrate that the Spitzer c2d surveys are complete 
to young objects with luminosities as low as O.O5L0. 

Finally, a comparison of our source list to the ded- 
icated search for ve ry low luminosity protostars by 
iDunham et al.l (|2008l ) confirms that we are not missing 
any embedded protostellar sources in Perseus down to 
the completeness limits of that survey {L^oi ^ 0.03 — 
0.05Lq). 

4. SOURCE CLASSIFICATION 

To study the early evolution of protostars, it is nec- 
essary to identify the evolutionary state of the embed- 
ded candidates in our sample. This is typically ac- 
complished by classifying sources into discrete groups 
based on SED characteristics. A number of classifi- 
cation methods are employed in the literature; most 
often used are the near- to mid-infrared spectral in- 
dex aiR = dlog(AFA)/d log(A) ( Lada 1987), the bolo- 
metric temperature T^oi (jMvers fc Laddil993.> see 
and the ratio of submillimeter to bolometric luminos- 
ity Lsubmm/Lboi ( Andre_eLinil993l) . Generally, Lsubmm, 
is taken to be the integrated luminosity at wavelengths 
A > 350 /im. 
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Fig. 5. — Distribution of Tjoj as a function of the distance to the nearest 1.1 mm core, for candidate embedded protostars in all three 
clouds. Distances are in units of the 1.1 mm core FWHM size (9imm)i and symbols are as in Figure |4] Sources within 1.0 X 6imm of a 
1.1 mm core position are considered embedded within that core. There is a clear correlation between smaller distances, or more embedded 
sources, and lower T;,;,; values, suggesting that Tf,o; is a good measure of evolutionary state for deeply embedded sources. For reference, 
the total numbers of sources in each quadrant of the plot are given in parentheses. 
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Fig. 6. — Distribution of am as a function of the distance to the nearest 1.1 mm core for embedded protostar candidates in all three 
clouds, similar to Figure[5] While there is some correlation between higher am values and smaller distances, it is not as compelling as the 
correlation observed for T^oj (Figure [Sjl. For reference, the total numbers of sources in each quadrant are given in parentheses. 
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Any of these methods must come with the caveat that 
protostellar mass may affect the cl assification, which onl y 
detailed modeling can resolve fe.g. lHatchell et al.ll2007aD . 
We typically have only one flux measurement for A > 
160 /im, so an accurate determination of Lsubmm is not 
feasible with these data.^ We will focus, therefore, on 
Tboi and aiR. 

4.1. Comparing Classification Methods 

In Figure |4] we compare classifications based on am 
and Thoi for the candidate embedded protostars in 
Perseus, Serpens, and Ophiuchus. Thick symbols indi- 
cate sources associated with 1.1 mm emission, while thin 
symbols denote sources with upper limits at 1.1 mm. The 
spectral index aju is determined from a least squares fit 
to all detections between A = 2/im and A = 24 /im. 
The calculation of Ttoi is descri bed in ^3.11 and We 
adopt Ttoi class divisions from IChen et a l. (199 5J) (see 
Table [T]) , and aiR divisions from lAndre fc Montmerld 
([1994): aiR < (Class I) and am > (Class II). No 
well-defined am criteria exists for Class sources, as 
deeply embedded objects were generally not visible in 
the mid-IR prior to Spitzer. 

The two classification methods agree fairly well for 
Class II and "warmer" Class I sources. Sources with 
Tboi ^ 100 K, by contrast, have a wide range of am 
values, and a few of these "cold" sources even fall into 
Class II based on the spectral index. One of the rea- 
sons for this large scatter in am becomes apparent when 
examining the SEDs of some deeply embedded sources, 
which show considerable differences at short wavelengths. 
Many protostellar SEDs are not monotonically increas- 
ing from 3.6 to 24 /zm, often falling from 5.8 to 8 /zm and 
rising again at longer wavelengths (e.g., Per-emb 18, Per- 
emb 22, Oph-cmb 1; Figured]). Geometric effects such as 
scattered light from an outfiow cavity or absorption from 
ices in the protostellar envelope are the most likely cause 
of these features. In any case, non-monotonic behavior 
at short wavelengths will clearly bias the calculation of 
am- Calculating am from a straight line fit between 2 
and 24 ^m, rather than a least squares fit, still results in 
a large range of am, values at low Tboi. 

To determine whether am or T{,oi is a more accurate 
measure of the true evolutionary state, we look at the cor- 
relation between both measures and the degree to which 
a given source is embedded. Figure O shows the distri- 
bution of Tboi with respect to the distance from each 
protostellar candidate to the nearest 1.1 mm core, in 
units of the core FWHM size. In all three clouds, es- 
sentially all sources with Tboi ^ 200 K are located within 
one core FWHM of a 1.1 mm peak. Thus the majority 
of "cold" objects, as defined by low Tboi, appear to be 
embedded in dense envelopes, and are likely to be at an 
early evolutionary stage. Furthermore, there is a clear 
correlation between smaller distance and lower Tboi , sug- 
gesting that sources with lower Tboi sue more embedded 
than those with higher Tboi- This correlation remains 
even if we exclude the 1.1 mm fiux from the calculation 
of Tboi, eliminating the possible bias between millimeter 

^ It is possible to calculate Lsubmm by assuming a modified 
blackbody spectrum and a value for /3 (e.g. Hatchell et al. 2007a), 
but given tfie assumptions involved we choose not to pursue this 
method. 



cores and Tboi- Thus, if we affiliate the degree of embed- 
dedness with youth, then Tboi appears to correlate well 
with evolutionary state. 

Figure [6] is similar, for am', while there is some corre- 
lation between higher am. values and smaller distances 
to the nearest 1.1 mm core, the relationship is not nearly 
as clear as for Tboi- In particular, there are a number of 
sources with large am values that do not seem to be em- 
bedded (the distance to the nearest core is much larger 
than one FWHM), and there are several deeply embed- 
ded sources with low ajR values. Based on Figures HHS] 
and a visual examination of sources with Tboi < 70 K, 
a number of which are known Class sources, we con- 
clude that Tboi is more reliable than am for classify- 
ing deeply embedded protostars. Hereafter, Tboi will be 
used to charac terize protos t ellar s ources, and we use the 
divisions from IChen et all (|1995D and Table [T] to place 
sources into Class 0, Class I, or Class II. 

5. CHARACTERISTICS OF CLASS 0, I, AND II SOURCES 

The total numbers of Class and Class I sources in 
each cloud, as defined by Tboi classifications, are given in 
Table m Note that in addition to Tboi, we also require 
Class I sources to be detected at 1.1 mm. This is not 
an unreasonable requirement, as our definition of Stage I 
dictates that they should be detected at millimeter wave- 
lengths {Menv > 0.1 Mq). The individual spectra (for 
A = 1.25 — 1100 /im) of all Class and Class I sources 
are shown in Figure [T] 

Numbers in parentheses in Table [5] are based on clas- 
sifications using the prismoidal, rather than midpoint, 
method for calculating Tboi; these give some idea of 
the uncertainties in the number of Class and Class I 
sources. Despite the sometimes substantial difference in 
the midpoint and prismoidal Tboi values for individual 
sources, there is very little difference in the resulting 
number of Class and Class I protostars, with the ex- 
ception of Serpens where the number of Class sources 
increases from 10 to 14. 

Statistics for Class II are not given; as we intended to 
select against sources without a protostellar envelope, we 
are necessarily incomplete to these objects. For example, 
there are 46 sources in our origi nal Serpens candidate 
sample with Tboi > 650 K, while iHarvev et all (|2007aD 
find 132 Class II sources. We do detect a few Class II 
sources at 1.1 mm, however (see i j5.ip . 

The average spectra of Class 0, Class I, and Class II 
sources in Perseus and Serpens are shown in Figure [51 
To calculate the average spectra, individual SEDs are 
weighted by 1/Lboi, so that we are not biased by the 
most luminous sources. Error bars in Figure [5] repre- 
sent the la error in the mean {crx/^/N\), but the dis- 
persion in the sample (cta) is much larger. Source SED 
shapes do not fall into discrete bins, but form a con- 
tinuous distribution between the averages shown. Error 
bars are large at wavelengths where many sources are 
not detected (at short wavelengths for Class 0, long wave- 
lengths for Class II, and 160/zm for all classes), and where 
there are significant variations from source to source. 

If Ophiuchus sources are included in the average calcu- 
lation, the resulting average spectra are skewed toward 
having more flux at short wavelengths, have larger dis- 
persion, and the similarity of SEDs within each bin is 
reduced. This behavior suggests that Tboi may be biased 
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Fig. 7. — Spectral energy distributions of all embedded protostars in Perseus, Serpens, and Ophiuchus, including photometry from 
2MASS (1.25, 1.65, 2.17 ^tm), IRAC (3.6, 4.5, 5.8, 8.0 Mm ), MIPS (24, 70, 160 Mm), an d Bolocam (1.1 mm) data. Also included are 
SCUBA 450 and 850 Mm fluxes from the literature ([Sandell fc Kneal200ll: IKirk et al.ll200g) and IRAS 25 or 60 Mm fluxes where they are 
used to replace saturated 24 or 70 Mm fluxes (see § 13.11 1. A sample is shown here; the full figure is available online only. 
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Fig. 8. — Average spectra of sources in Class 0, Class I, and Class II, for the Perseus and Serpens samples. Sources are classified based 
on bolometric temperature, with the additional requirement that Class I sources be detected at 1.1 mm (see Table[TJ. Individual SEDs are 
weighted by l/Li,oi in the average calculation. Error bars are the Icr error in mean, and not the sample dispersion, which is much larger. 
The Class II spectrum is not necessarily representative of all Class II objects, as we are very incomplete to sources without protostellar 
envelopes. 



for many sources in Ophiuchus, hkely due to the lack of 
information at 160 /im, so we exclude Ophiuchus from 
the average spectra. 

The progression from Class to Class I to Class II 
is consistent with a sequence of physical evolution. In 
particular, the average Class spectrum has the low- 
est flux densities from 1.25 — 24 /im, as expected for 
deeply embedded sources with massive, extincting en- 
velopes, and the highest fluxes at 70 — 1100 /Ltm, where 
reprocessed protostellar flux is emitted by the cold en- 
velope. By contrast, the average Class II spectrum is 
relatively flat, with a much larger percentage of the pro- 
tostar flux emerging at shorter wavelengths, as expected 
for older sources without much circumstellar material. 
One must keep in mind, however, that Ttioi is defined 
such that spectra which peak at longer wavelengths will 
have a lower Tboi. 

Large error bars on the average Class spectrum for 
A = 1— 3 fim. are indicative of the widely varying behavior 
of Class objects in the near-infrared, and the presence 
of non-monotonic behavior is apparent at A = 3.6 ^m. 
The large separation between the Class and Class I 
spectra, as well as the continuous range of SED shapes, 
suggests that a more continuous means of estimating evo- 
lutionary status is preferable to the standard classes. 



5.1. Class II Sources with 1 mm Emission 

The average Class II spectrum is not necessarily rep- 
resentative of all Class II objects, due to severe incom- 
pleteness to non-embedded sources. The Class II sources 
we do detect are likely relatively young Stage II objects, 
before a substantial fraction of disk mass is dispersed or 
accreted. 

Given our 1.1 mm sensitivity limit of approximately 
0.1 M0, almost all sources detected at 1.1 mm will be 
dense cores or envelopes around relatively young proto- 
stars. In general, by the time a protostar has consumed 
or dispersed its massive e nvelope and enters S tage II (or 
the T Tauri phase, e.g., I Adams et al.lfl987f ). the rem- 
nant disk of gas and dust has too little mass to be de - 
tected by our millimeter surveys. iCrapsi et all (l2008h 
define Stage II as sources with a circumstellar mass be- 
low 0.1 M0. Typical me asured masses of Class II disks 
are 0.01 - 0.1 Mq (e.g., iBeckwith fc SargentI [19961 ) . al- 
though values as larg e as 1 have been measured 
(jBeckwith et al.lll990f ). 

In a few a cases we do detect 1.1 mm emission around 
sources with Class Il-type SEDs. These objects have 
bolometric temperatures T^oi > 650 K and a flux density 
at 1.1 mm that is lower than the flux densities from 3.6 
to 24 /im. With the exception of two "flat spectrum" ob- 



12 



Enoch et al. 



jects (-0.3 < aiR < 0.3: lGreene et al.lll994f) . the near- 
to mid-infrared spectral indices, aiR, of these objects are 
in the range —0.34 to —1.04, confirming their Class II sta- 
tus. In some cases, the 1.1 mm emission is unresolved, 
consistent with a compact disk. Often, however, these 
Class II sources are in a region of confused millimeter 
emission, so their physical association with the 1.1 mm 
emission is not secure. 

Table [U lists all Class II sources (650 < Tboi < 2800 K) 
in each cloud that are detected in our 1.1 mm surveys. 
Those in confused regions of 1.1 mm emission are indi- 
cated by a "*" in the "Bolocam ID" column. Sources 
with point-like 1.1 mm emission emission centered on 
the Spitzer position (3 in Perseus, 1 in Serpens, and 3 in 
Ophiuchus) may have massive disks. 

5.2. Individual Sources 

Here we briefly discuss a few examples of newly- 
identified or otherwise interesting embedded protostellar 
sources. 

5.2.1. IRAS 03292+3039 

IRAS 03292-^3039 (Per-emb 2) is a little-studied. 
Class source associated with the 1.1 mm core Per- 
Bolo 66; it was discussed briefly in Paper I. An image 
of IRAS 03 292+3039, together wit h the SED, is shown 
in Figure m iJorgensen et al.l ()2006l ) identified this as an 
outflow source, noting the large-scale outflow visible in 
the 4.5 /im IRAC band. 

Bright 1.1 mm emission centered on this object indi- 
cates a fairly massive protost ellar envelope (2.9 M(7i) , 
while the powerful outflow (IWalawender et al.l 120051 : 
iJargensen etall l2006t iHatchell et al.l l2007bl ) and low 
bolometric temperature {Ttoi = 25 K) are evidence 
of an extremely young, energetic embedded protostar. 
The spectral energy distribution is similar to well-known 
Class protostars in Perseus such as NGC 1333-IRAS 4, 
and nearby IRAS 03282-^3035. The smaU fan-shaped 
nebulosity visible at 3.6 and 4.5 /zm (Figure [9]) is most 
easily explained by a cone-shaped cavity, carved out of 
the dense envelope by an energetic outflow. The orien- 
tation of the nebulosity corresponds w ell to the larger- 
scale outflow traced by IRAC emission (jJorgensen et al.l 
l2006l ) , which is indicated by thick gray lines in Figure M 
The one-sided nebulosity at 3.6 /xm is strikingly similar 
to the Whitney et al. (2003) model IRAC image of an 
early Stage sources viewed at an inclination angle of 
30° (their Figure 12a). 

5.2.2. Per-Bolo 102 

Per-Bolo 102 is a bright 1.1 mm source t hat was iden- 
tified in the Bolocam survey of Perseus ([Enoch et al.l 
I2006D . It lies within the region of active star-formation 
near IC 348, w hich includes the fanious o utflow-driving 
source HH 211 (jMcCaughrean et al.lll994f ). The Spitzer 
24 /xm image resolves the luminous internal source into a 
double object (Per-emb 16 and Per-emb 28). While the 
two SEDs are very similar for A > 8 /im (Figure [9]), vari- 
ations at shorter wavelengths cause Per-emb 16 (T&oj = 
56 K) to faU into Class 0, while Per-emb 28 {Uoi = 72 K) 
falls just outside the Class 0/Class I boundary. These 
sources are good examples of why a more continuous 
evolutionary scheme is preferable to the standard clas- 
siflcations. 



Recently obtained CARMA interferometric observa- 
tions at A — 2.7 mm resolve the millimeter core into 
two sources, coincident with the Spitzer source positions 
(M. L. Enoch et al., in preparation). The flux ratio at 
2.7 mm of the northern (Per-emb 16) to southern (Per- 
emb 28) source is at least 2:1, further evidence that Per- 
emb 28 is slightly more evolved. Although not as massive 
or cold as nearby HH 211 (Per-emb 1) and IC 348-mms 
(Per-emb 11), Per-Bolo 102 is an interesting case study. 
It may be a separate-envelope binary system, with two 
nearly coeval Class or early Class I sources. The sources 
are separated by 17", or 4200 AU. Binary separations of 
this order are consistent with early fragmentation in a 
relat ively dense cloud ("prompt initi al fragmentation" , 
e.g.. |Pringldll989HLoonev et al.ll2000f) . in which case the 
individual sources would have distinct protostellar en- 
velopes. 

In a binary formed via gravitational fragmentation, we 
would expect the separat ion to correspond to the local 
Jeans length ()Jeanslll928h : 

/ TTC^ \ 

where Cs is the local sound speed, and fip = 2.33 and n 
are the mean molecular weight and mean particle density, 
respectively. A Jeans length of 4200 AU would require a 
relatively high density (n ~ 6 x 10^ cm^'^, assuming Cs = 
0.2 km s~^). The mean density of the Per-Bolo 102 core, 
measured within an aperture of 10^ AU, is 4 x 10^ cm~'^, 
close to the required value. 

5.2.3. Ser-Bolo 33 

One noteworthy example of a Class II source 
associated with 1 mm emission is Ser-Bolo 33 
(SSTc2dJ183006.12-f 004233.8) in Serpens (Figure HD, a 
very bright Spitzer source with Lhoi — 3.6 Lq, Tf,oi = 
871 K, and compact 1.1 mm emission centered on the 
Spitzer position. The spect ral index, arR = ~ 0-42, also 
places this object in Class Il. lVieira et al.l(|2003[ ) included 
this source in a sample of Herbig Ae/Be candidates, be- 
lieved to be the intermediate mass Class II counterparts 
of low-mass T Tauri objects, although the spectral type 
is F3. The 1.1 mm mass calculated assuming an optically 
thin disk is 0.17 M©, approximately 10% of the stellar 
mass, ~ 1.3 M©, estimated fr om the measured effec- 
tive temperature (Tos ~ 6300 K. IVieira et al. 2003) and 
an empirical Toff — M relation (|Habets fc Heintze.1981. ). 

6. ALTERNATIVE CLASSIFICATIONS 

One advantage of such a large sample is that it allows 
us to define a more continuous evolutionary sequence 
than the standard classes that were appropriate for the 
smaller samples previously available. With this in mind, 
we divide the protostellar sources in each cloud into 
smaller Tj^oi bins than those of the standard Class O/I/II 
divisions. Average spectra for "early Class 0", "late 
Class 0", "early Class I", "late Class I", and Class II 
(see Tabled] for definitions) sources in Perseus and Ser- 
pens are shown in Figure [Till Error bars represent the la 
error in the mean, and average spectra are calculated as 
described in the previous section. The largest error bars 
are seen for the shortest wavelengths in the lowest Tboi 
bins, where NIR fluxes vary significantly from source to 




Fig. 9. — Three-color Spitzer images and SEDs of individual sources discussed in i|5.2l Images are (8,24,70 fim), unless otherwise 
noted, with 1.1 mm contours. Left: Three-color (3.6,24,70 fim) image of IRAS 03292+3039 (Per-emb 2). Fan-shaped nebulosity is visible 
at 3.6 ^m, most likely scattered light from a narrow outflow cavity in this deeply embedded Class source. Thick lines indicate the 
approximate orientation of the larger-scale outflow (j0rg ensen et ai. .! 2d0(j ). Center: Per-Bolo 102, a new candidate binary Class O/I source. 
Although the SEDs of the two embedded protostars look very similar at A > 8 fim, the southern source is brighter at shorter wavelengths, 
fainter at 2.7 mm (see text), and has a slightly higher Tjoj. Right: Ser-Bolo 33, a bright Class II object (Table|6]l. While this is clearly a 
more evolved source, it is massive enough that we detect 1.1 mm emission from a compact disk. 
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Fig. 10. — Average spectra in designated Tjoj bins (see Table[T]l, using sources from Perseus and Serpens. The average is calculated as 
described in Figure |8] and the text, and error bars represent the error in the mean. The progression in SED shape for increasing T^oj is as 
expected if this represents physical evolutionary sequence. 



source, likely depending on outflow opening angle and 
viewing geometry. Binning of Tboi is based on general 
agreement of SEDs in a given bin, as determined by eye, 
and should not be interpreted as strict boundaries. 

While there are significant similarities between the 
SEDs within each bin, individual source SEDS fill the 
continuum between the average spectra, and each aver- 
age spectra lies within the la dispersion of neighboring 
Tboi bins. The average spectra change as one would ex- 
pect if an extincting envelope is gradually accreted or dis- 
persed, with the protostar becoming more visible at short 
wavelengths. At all wavelengths except 24 and 70 /im, 
the flux rises or falls monotonically with increasing T^oi . 
For these intermediate wavelengths, the observed flux 
may rise initially as hotter dust close to the protostar 
is revealed, then fall as the mass of circumstellar mate- 
rial drops. Despite the relatively narrow bins, there is 
still a rather large change between "late Class 0" and 
"early Class I", particularly at A = 3.6 — 24 /im; this 
transition may occur rapidly, or these wavelengths may 
be particularly sensitive to geometry. 

7. COMPARISON TO MODELS 

Even for inflnitely well-sampled SEDs, spectrum shape 
is not necessarily directly correlated with age or de- 
gree of embeddedness. For example, viewing geometry 



can have a strong effect on SED shape . In three di- 
mensional rad iative transfer models of (IWhitnev et al.l 
,2003,: Robitaille et al.ll2006t ICrapsi et all 120081 ). Stagel 
and Stage II sources can have quite similar SEDs when 
viewed at the right inclination angle (e.g., when the ob- 
server's line-of-sight intersects the outflow opening angle 
of Stage I sources). Even Stage sources can appear 
much warmer if we happen to be looking directly into 
the outflow cavity, although the probability of that oc- 
curring is small. 

We compare our spectra to the results of protostel- 
lar models, which predict protostellar spectra based on 
source age, mass, accretion rate, etc., both to gain in- 
sight into the evolutionary state of sources and to eval- 
uate the effectiveness of such models in matching ob- 
served sources. Rather than model each source individu- 
ally, we compare the average spe ctra fro m Figure fTUl with 
predicted spectra from Whitnev et al.l ([2003,) for "early 
S tage 0" , "late S t age 0" sources, etc. in Figure fTTl 

IWhitnev et al.l (|2003D begin with a set geometry for 
each stage, including some combination of accreting pro- 
tostar, flared protostellar disk, infalling envelope, bipolar 
outflow, and grain models for each region, then use ra- 
diative transfer modeling to predict protostellar spectra. 
Colors in Figure [11] correspond to inclination angles of 



Protostars in Serpens, Perseus, and Ophiuchus 



15 



XF,^ (ergs s ^ cm ^) 
o o o o 




o 

^ o 
o 



o 
o 
o 



XF^ (ergs s ^ cm ^) 
o o o o 









Q 
















\^ r-t- 























o 



AF^ (ergs s ^ cm ^) 
o o o o 




o 



XFjj (ergs s ^ cm ^) 
o o o o 





Fig. 11. — C omparison of the average observed spectra from Figure [TPl with predicted spectra from the protostellar evolution models of 
IWhitnev et al.l {2003) (plotted as uSi/). Distinct geometries are chosen to represent five evolutionary stages, and include a combination 
of an accreting protostar, a protostellar disk, envelope, and a bipolar outflow. Colors indicate inclinations angles of the protostellar disk, 
from 0° (pink) to 90° (dark green). With some exceptions the Stage and Stage II models match the observed spectra fairly well, but the 
Stage I models over-predict the near-IR and under-predict the mid-IR flux. In general the models tend to underestimate the millimeter 
flux and thus the mass of the protostellar envelope. 
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the protosteUar disk from 0° (pink) to 90° (green) , where 
the outflow is perpendicular to the disk. Model envelope 
infall rates decline from 10~^ MQyr"'^ in early Stage to 
10~® M0yr~^ in late Stage I. Similarly, the disk radius 
and cavity opening angle increase and the cavity density 
decreases as one moves from early Stage to late Stage I. 
Both models and average observed spectra (thick black 
lines) are scaled to a total luminosity of ILq. 

The model early Stage and observed early Class 
spectra agree fairly well. The observed spectrum lies 
above most of the models at A < 24/im, which could 
be explained by gaps in the inner envelope of some 
protostars, allowing s hort wavelength flux to escape 
(jJorgensen et al.l l2005f). or a wider outflow angle than 
the model (5°). The IWhitnev et al.' (2003) models pre- 
dict that a small fraction (approximately 1/10) of early 
Stage sources will have spectra similar to Stage I, at 
very low inclination angles (looking down the outflow). 
The observed Class II spectrum is also consistent with at 
least the low inclination models; as we select for sources 
with 1.1 mm emission it is not unexpected that the ob- 
served millimeter point is higher than predicted by the 
models. For late Stage the agreement is again pretty 
good, except that the model under-predicts the millime- 
ter flux, and thus the mass of the protosteUar envelope, 
a feature present in all spectra later than early Stage 0. 

By contrast, the agreement between both Stage I mod- 
els and the observed average spectra is quite poor. The 
models severely over-predict the near-IR flux and under- 
predict the mid-IR flux. In fact, the characteristic 
double-peaked profile of the Stage I models is almost 
n ever seen in the ob s erved spectra (Figure [T]). 

Ilndebetouw et al.l ()2006l ) demonstrate that envelopes 
with clumpy density distributions can increase the ob- 
served flux at 10 /im, with short wavelength emis- 
sion escaping through low opacity regions between the 
clumps. This kind of clumpy distribution eliminates 
some of the "double- peaked" structure apparent in the 
IWhitnev et"an (|2003[ ) models, which assume a smooth 
density profile in the envelope and constant density in 
the outflow cavity, although the clumpy models also in- 
crease the flux at < 3/im. Alternatively, the presence of a 
larger, thicker protosteUar disk could flatten the Stage I 
spectra by absorbing near-IR flux and producing more 
mid-IR flux. A foreground cloud could also absorb near- 
IR fiux, without producing any mid-IR flux. 

It seems that we still do not understand how sources 
transition from Stage to Stage I, at least in relation 
to these models. Whether it is the physical model that 
is unrealistic or how that physical geometry translates 
to the observed spectra is not clear. It is important to 
note that the same a uthors have deve loped a more com- 
plete grid of models (jRobitaille et al.i ,2006): we compare 
here to the earlier models rather than fitting each av- 
erage spectrum because we want an idea of the global 
agreement with evolutionary stage. 

8. PROTOSTELLAR EVOLUTION 

8.1. Luminosity Evolution 

We examine the evolution of embedded protostars in 
more detail using the Ltoi — Thoi diagram, the proto steUar 
equivalent of the H-R diagram (" Mvers et al.l [l998^. Here 
Tfjoi is used as a measure of temperature rather than T^ff, 
which is not well-characterized for embedded sources. 



As discussed in I M vers et al.l (|1998D , newly-formed pro- 
tostars should begin at low Lf,oi and Ti,oi, increasing in 
both Lboi and T^oi as accretion proceeds. If accretion 
abates or is otherwise halted, then L^oi will decrease for 
steadily increasing Tboi- Eventually, sources will move 
onto the main sequence, at Ti,oi > 3000 K. 

Figure [T^j plots Ltoi versus Tboi for the candidate em- 
bedded protostar samples in Perseus, Serpens, and Ophi- 
uchus. Filled and open symbols are used for sources that 
are associated with 1.1 mm emission, while sym- 
bols indicate sources with upper limits at 1.1 mm.^ All 
sources detected at 1.1 mm are further divided accord- 
ing to whether they have reliable measured 160 /xm fluxes 
(squares), are saturated at 160 ^m (circles), or are not 
detected at 160 /im (triangles). Given the discussion in 
€A1 regarding the effect of missing 160 /im fluxes on the 
calculation of Tboi, circles would be expected to move 
up and to the right in this diagram for "cold" sources 
{Tboi ^ 100 K), and down and to the left for "warm" 
sources {Tboi ^ 100 K), if 160 /xm fluxes were available. 

While upper limits at 1.1 mm ("-I-" symbols) are ex- 
pected for more evolved sources with Tboi ^ 500 — 600 K, 
colder sources with no 1.1 mm detection may be ei- 
ther misclassifled (e.g., background galaxies), or very 
low mass sources whose 1.1 mm flux is below our de- 
tection limit. Given our physical definition of Stage I 
{M* ^ Menv', Menv > O.IM©), embedded sources not 
detected at 1.1 mm must have stellar masses less than a 
few tenths of a solar mass, or have very little remaining 
envelope (M* ^ Menv)- 

ProtosteUar evolutionary tracks, which predict source 
properties as a function of age and mass, can easily be 
compared to our data using t h e Lbo i — Tboi diagram. 
Model tracks fro m iMvers et all (Il998fl ( solid line; here- 
after M98) and lYoung fc EvansI (| 2005l) (dotted lines; 
hereafter YE05) are shown in Figure [121 YE05 adop t 
the standard inside-out collapse model of S hu! (|1977t l. 
and assume that no mass is lost in the formation pro- 
cess. From top to bottom, the YE05 models are for 
sources with masses of 3.0, 1.0, and 0.3 Mq, assuming 
a constant accretion rate o f dM/dt = c^/G, where Cg is 
the effective sound speed (|Shulll977l ). From the initial 
singular isothermal sphere, finite masses are achieved by 
truncating the outer radius of the envelope. A one di- 
mensional radiative transfer model (DUSTY) is used to 
calculate observational signatures {Lboi, Tboi, etc.) from 
the accretion model. 

Unlike YE05, M98 do not assume that that the entire 
mass of the original core ends up in the final star, but 
rather that a significant fraction of the core mass is lost 
in the star formation process. The M98 model shown is 
for a source with initial core mass of 1.8 Mq and final 
stellar mass of 0.3 Mq. M98 assume an accretion rate 
that is initially dM/dt = C^/G, but falls off exponentially 
with time, designed to match the observed luminosity of 
pre-main sequence stars at Tboi ^ 3000 K. Thus, the lu- 
minosity is significantly lower than the YE05 tracks at 
later times. Both evolutionary models assume an accret- 
ing central protostar, a circumstellar accretion disk, an 
extended envelope, and a contribution to the luminosity 

^ Note that although we plot sources with 1.1 mm upper limits 
here, our final embedded protostar samples (Tables |2]jiJ includes 
only sources with 1.1 mm detections. 
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Fig. 12. — Bolometric luminosity versus bolometric temperature (Lfcoj — Tfj^i) diagram for embedded protostar candidates in Perseus, 
Serpens, and Ophiuchus. Filled and open symbols indicate that a given source is associated with 1.1 mm emission, while "+" symbols 
indicate upper limits at 1.1 mm. Symbol are further divided based on whether there is a reliable 160/^m measurement, an upper limit, or 
if we are unable to measure a 160/^m flux due to saturation. The right axis shows the mass accretion rate, dM/dt, calculated from L^oj 
using Eq. |6]and assuming Af» = 0.5Mq and i?* = 5R p,. Model evo lutionary tracks from Young & Evans (2005) for sources with 3.0, 1.0, 
0.3 Mq (dotted lines, from top to bottom) and from lMvers et all |(T99^ ) for a source with final mass 0.3 Mq (solid line) are shown for 
comparison. The large population of low-luminosity (Lf,oi < O.IL0) Class I objects argues strongly for episodic accretion by Stage I. 



from gravitational contraction of the protostar. YE05 
also include nuclear (Deuterium) burning. 

Ignoring for a moment the population of Class I sources 
in each cloud with Lhoi values well below both mod- 
els, the M98 model, for which a large fraction of the 
core mass is ejected or dispersed, is more consistent with 
the observed protostellar sources. The ratio of the fi- 
nal stellar mass to initial core mass for the M98 model 
shown (/eff = Mcore/Mstar = 0.3/1.8 = 0.1 7) is Smaller 
than the v alues (/sfF = 0.3 ib .1) found bv lAlves et al.l 
(|2007f) and lEnoch et cil] (|2008f) by comparing the shape 
of the core mass distribution to the initial mass function 
(/eff = 0.3 ± 0.1 and fed ^ 0.25, respectively). 

The slightly better match to the M98 model may be 
irrelevant, however, as neither a constant nor a steadily 
decreasing accretion rate is consistent with the observed 
protostellar populations in these clouds; many sources lie 
below all four model tracks. Our data con firm and exac- 
erbate the "luminosity problem" noted by iKenvon et al.l 
(fl99nh . that Class I protostars in Taurus were observed 
to have lower Li,oi values than expected based on the av- 
erage mass accretion rate required to make a IM0 star. 

In particular, the large population of Class I sources 
with low Lboi in each cloud is difficult to understand 
in relation to most existing protostellar evolutionary 
models. A general feature of such models is that the 
bolometric luminosity peaks in the Class I stage, a re- 
sult that is true for constant accretion rates (YE05), 
decreasing rates (M9 8), and gravo-turbulent models 
(jFroebrich et al.l[2006f l. In contrast, we find quite similar 
mean luminosities for the Class and Class I samples: 
2.4 L0and 2.2 L©, respectively. The standard devia- 
tion for both samples is large: 3.5 and 4.7 for 



Class and Class I, respectively. The median luminosity 
of Class I sources is a factor of 3 lower than the mean, 
0.7 Lq. 

8.2. Episodic Accretion 

Given the above discussion and the large observed 
spread in L^oi of three orders of magnitude for Class I 
sources, much larger than the range in envelope masses, 
we conclude that mass accretion during the Class I stage 
is episodic. Class I sources with low Ljjoi can be ex- 
plained by periods of relative quiescence when the bolo- 
metric luminosity, which is driven primarily by accre- 
tion luminosity, drops by at least a factor of 10. Con- 
versely, Class I sources with high Lboi values, e.g., those 
that form the upper envelope of the distribution and ap- 
pear to be consistent with the YE05 models, would cor- 
respond to periods where the accretion rate is close to 
the "Shu accretion" value. Approximately 20% of the 
Class I sources have L^oi < O.IL© (7/39 in Perseus, 5/25 
in Serpens, and 5/25 in Ophiuchus), and could be con- 
sidered candidate very low luminosity objects (VeLLOs; 
iDunham et"alll2008l: iDi Francesco et al.ll2007D . 

Episodic accretion is not an unreasonable solution; 
evidence for variable mass accretion and ejection is 
plentiful, includin g that based on mod eling FU Orio- 
nis eruptions (e.g.. Hartmann fc Kenv on 1985) , and bow 
shock s in Herbig-Haro outflows (e.g., Reipurth fc Ballvl 
|2001| ). In traditional models of episodic accretion, infall 
from the envelope onto the disk is constant and accre- 
tion from the disk onto th e protostar is episodic (e.g., 
iKenyo n fc IIart mannlll9 95l. for example due to g ravita- 
tional instabilities in the disk (Vorobvov fc Basul 12006.) . 
Other scenarios such as the "spasmodic" infall model of 
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Fig. 13. — Average spectra of low-luminosity Class I sources and "non-envelope Class I" sources (those that have Tjoj < 650 K but no 
1.1 mm emission). Average Class 0, I, and II spectra from Figure[8]are shown (dashed lines) for reference. The Iow-Ljqj Class I sources have 
lower short wavelength points relative to the 160 fim and 1.1 mm fluxes, as expected for low accretion rates, but are otherwise consistent 
with the average Class I spectrum. The average "non-envelope Class I" spectrum is intermediate between the average Class I and Class II 
spectra, as expected if these sources represent an transitional stage in which there is little or no remaining envelope. 



iTassis &: MouschoviasI (|2005l ). where material is held up 
by the magnetic field at the inner edge of the envelope 
are also plausible. 

Although the observed distribution of L(,oi for Class I 
sources argues strongly for episodic accretion, that is not 
the only possibility. One alternative is that the \ow-Li,oi 
Class I sources are simply very low mass objects, and 
we are somehow missing their low-mass Class counter- 
parts. Based on our sensitivity limits at 70, 160 /im and 
1.1 mm, we should be able to detect Class sources with 
Lboi ^ O.O5L0. The lower limit to the observed Class 
luminosity, however, may not be the internal luminosity, 
but heating of the envelope by the interstellar radiation 
field (ISRF), which can contribute as much as 0.2 — 0.3Lq 
to Lboi for envelope masses of 1 — 3Mq (YE05). For the 
same reason, a Class source with a very low mass accre- 
tion rate would not necessarily have a very low bolomet- 
ric luminosity, making it difficult to determine if episodic 
accretion is already present at the Class stage. 

In Figure I13[ we show the average spectrum of low- 
luminosity {Lfjoi < O.2L0) Class I sources in Perseus 
and Serpens. On average, these sources have suppressed 
short wavelength points relative the 160 /zm and 1.1 mm 
fluxes, but are otherwise consistent in shape with the 
average Class I spectrum.^ The average low-L;,,,; spec- 
trum is consistent with low accretion luminosities (evi- 

Note that the average spectra are normalized by l/Ljoj, so 
on an absolute scale the low-L^oj spectrum would be substantially 
fainter than the Class I spectrum. 



dent at short wavelengths) in sources with "normal" en- 
velope masses, as expected if the low luminosity objects 
are similar to other Class I sources but with lower mass 
accretion rates. The characteristics of low luminosity 
protost ars in Perseus, Serpen s, and Ophiuchus are ana- 
lyzed in lDunham et all (|2008D . 

8.3. Mass Accretion Rates 

Assuming that the bolometric luminosity in Class 
and Class I is due entirely to accretion, we can estimate 
the accretion rate from Ltoi ■ 



M 



dM 
~dt 



(6) 



where i?* and are the radius and mass of the em- 
bedded protostar. For the following we assume M» = 
O.5M0 and i?* = 5i?0, but accretion rates can be eas- 
ily scaled for different values of the mass or radius. 
A mass of O. SMfri is consist ent with the initial mass 
function (e.g IChabrieil I2003D and with ave rage YSO 
masses ()Merm et al.ll2008t ISpezzi et al.ll2008[ ): a radius 
of 3 — 5i?0 is typi cal for pre-main sequence models of low 
mass sources (cg lPallafc Stahler 1991: Robitail le et al.l 
l2006h . 

For comparison the "Shu accretion rate" of c^/G, 
which is used in the YE05 models, is approximately 
4 X 10-6 M0 yr-i for c, = 0.2 km s'^ The M98 model 
begins with M ^ 10^^ Mq yr~\ falling to 10"^ M0 
yr"^ by the time it reaches the main sequence. Making 
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a 1 Mq star in 5.4 x 10^ yr requires an average accretion 
rate of approximately 2 x 10~® Mq yr~^. 

The mean luminosity of Class I sources corresponds to 
M ^ 1 — 2 X 10^^ M0 yr~^, not far from the average 
required to make a solar mass star. By contrast, the 
luminosity of the low-Lboi {Lboi ~ O.IL©) Class I sources 
implies a mass accretion rate of only 7 x 10~^ Mq yr~^, 
substantially smaller than the Shu value or the average 
required to make a solar mass star. This result again 
suggests that either these sources will form very low mass 
stars (M* < O.OSM©)*, or that these sources are in a 
suppressed accretion phase of an episodic cycle. If we 
define "sub-Shu" accretion to be M < 10~^ Mq yr~^, or 
Lboi S ILq, then approximately 55% of Class I sources 
are observed to be in such a "sub-Shu" accretion state. 

If the episodic accretion picture is correct there must be 
a population of sources with accretion rates much higher 
than the average; indeed, the highest observed luminosi- 
ties {Lboi ~ 10 - 30 Lq) imply M ~ 1 - 2 x IQ-^ Mq 
yr~^, at least a factor of five higher than the Shu value 
or the average required to make a IMq star. Approxi- 
mately 5% of the Class I sources appear to be in such a 
"super-Shu" accretion rate {Lboi > IOLq, or AI > 10^^ 
Mq yr-1. 

Unfortunately, it is very difficult to determine the duty 
cycle of the episodic accretion without detailed star for- 
mation and accretion models. The small fraction of 
sources observed to have very high accretion rates (5% 
with M > 10~^ Mq yr~^), suggests t hat periods of rapid 
accretion must be fairly short lived. I Evans et al.l ()2008l ) 
employ a simple accretion model to estimate the time 
spent in accretion and quiescent phases, finding that half 
the mass of a 0.5Mo star could be accreted during 7% of 
the Class I lifetime. 

8.4. Envelope Mass Evolution 

Protostellar evolution models also predict the evolu- 
tion of envelope mass, Mf.nvi with T^oi, as shown in Fig- 
ure 1141 Symbols are similar to Figure I12[ with upper 
limits at 1.1 mm represented by arrows, and solid lines 
indicating the YE05 evolutionary tracks (3.0, 1.0, and 
0.3 Mq top to bottom). Some of the higher Tboi sources 
have high 1.1 mm upper limits because they lie in re- 
gions of extended or confused emission. Our determina- 
tion that these sources are not actually associated with 
the 1.1 mm emission is based on visual examination of 
the images and SEDs. 

Sources in all three clouds show a weak but consis- 
tent trend of decreasing Mf-nv with increasing T^oi , as 
expected if the envelope is gradually depleted by accre- 
tion onto the protostar. This trend is not a result of in- 
cluding the 1.1 mm flux, from which the envelope mass 
is derived, in the calculation of Tboi] the same trend is 
apparent even if the 1.1 mm point is excluded in the 
calculation of T^oi ■ 

Here, the YE05 model tracks fit the observed Menv — 
Tboi distribution quite well. Thus, a constant envelope 
infall rate reproduces the decrease in envelope mass with 
increasing Tboi, although it does not fit the evolution of 
Lboi with Tboi (Figure 

* Although this would require these sources to have a quite low 
core-to-star formation efficiency of < 10%, as only a few Class 
sources have envelope masses less than 0.5Mq . 



For standard episodic accretion models, where in- 
fall from the envelope is steady and the accretion rate 
from the disk to the pr otostar is the variable quantity 
(e.g.. lKenvon fc Hartma nn 1995), the envelope mass will 
steadily decrease with increasing Tboi even for variable 
accretion. Note that accretion stops when the enve- 
lope mass, as defined by the outer radius, has been ex- 
hausted. In the context of the YE05 models, the spread 
of Menv as a function of Tboi suggests stellar masses in 
the range M ~ 0.3 — 4 Mq in Perseus and Serpens, and 
M ~ 0.1 - 1.5 Mq in Ophiuchus. 

9. LIFETIME OF THE CLASS PHASE 

The length of time that sources spend in the Class 
phase is an important diagnostic of protostellar evolu- 
tion and how accretion rates evolve with time. If the 
rate of star formation in these clouds is steady in time 
(i.e., not occurring in bursts), and if Class 0-Class I- 
Class II represents a true evolutionary sequence, then 
we can use the number of objects in consecutive evolu- 
tionary phases to estimate the relative lifetimes of those 
phases: ti/t2 = Ni/N2. As we ultimately calibrate life- 
times based on the lifetime of Class II disks (~ 2 x 10^ yr ; 
iKenvon et aI1ll990t ICieza et al.ll2007l : ISpezzi et al1l2008f l , 
star formation must have been steady in time for at least 
the last 2 Myr. In addition, there can be no signifi- 
cant dependence of the evolutionary timescales on source 
mass. 

It is unlikely that all of these assumptions hold in ev- 
ery star forming region. In fact there is some observa- 
tional evid^n£e_Joi^mBSS-dependent evolution of dense 
cores (iHatchell et al.|[2008[) . which could easily translate 
into mass-dependent evolution after protostar formation. 
We can, however, hope to mitigate the effects of any 
breakdowns in our assumptions by utilizing our large 
sample and averaging over three different env ironments. 

In a companion paper ()Enoch et all l2008l ) we use a 
similar argument to derive the lifetime of the prestellar 
phase from the ratio of the number of starless cores to the 
total number of embedded protostars (Class -I- Class I). 
There we find there that the lifetime of dense {n > 10"* 
cm~^) prestellar cores is approximately equal to the life- 
time of the embedded protostellar phase, or 2 — 5 x 10^ yr, 
suggesting that such cores evolve dynamically over a few 
free-fall timescales. 

The relative number of Class and Class I sources in 
each cloud, and the Class lifetime derived from that 
ratio, is given in Table [5] Given in parentheses are the 
values resulting from using the prismoidal, rather than 
midpoint method for determining Tboi] the Class life- 
time is based on an average of the two methods, and the 
difference between them, added in quadrature with \/N 
statistical uncertainties, yields an uncertainty for the life- 
time. 

There are approximately half as many Class 
as Class I sources in both Perseus and Serpens 
(^ciass o/^ciass I = 0.7 and 0.4, respectively), suggest- 
ing that the Class phase lasts roughly half as long as 
the Class I phase. We adopt a total embedded phase 
hfetime of t^mb = tern + icii 5.4 x 10^ yr, derived 
based on the rela tive number of emb edded protostars and 
Class II sources vans "ei~an[2008h . Thus our measured 
ratios imply Class lifetimes of iciass o ~ 2.2 x 10^ yr in 
Perseus, and 1.7 x 10^ yr in Serpens. In Ophiuchus there 




Fig. 14. — Envelope mass versus bolometric temperature for embedded protostars in Perseus, Serpens, and Ophiuchus. Filled and open 
symbols indicate that a given protostar is associated with 1.1 mm emission, w hile sources with upper limits at 1.1 mm are plotted as 
arrows. Solid lines show the pr edictions of prot o stella r evolutionary models from lYoung fc EvansI |(S)05) for sources of mass (from top to 
bottom) 3.0, 1.0, 0.3 Mq. The lYoune &: EvansI II2005I) models describe the evolution of envelope mass with T^oj quite well, for sources of 
mass M ~ 0.3 — 4 Mq in Perseus and Serpens, and M ^ 0.1 — 1.5 Mq in Ophiuchus. 



are only 3 or 4 Class sources, resulting in a ratio of 
A^ciass o/A^ciass I = 0.1 - 0.2 and ^ciass 0.7 X 10^ yr. 

Taking all three clouds together, we find a lifetime for 
the Class phase of 1.72 ± 0.25 x 10^ yr. This value is 
significantly longer than a number of previous estimates 
of iciass ^ 10"* yr, based both on the number of Class 
sources in Ophiuchus (Andre & Montmcrlc 1994), and 
on comparison to evo lutionary models (2 — 6 x 10^ yr; 
iFroebrich et alJ 120061 ) . A short Class hfetime is gen- 
erally interpreted as evidence for a period of very rapid 
accretion early on in the evolution of protostars, caus- 
ing them to quickly reach Stage I, at which point the 
accretion rate decreases significantly. Our results argue 
against such a rapid accretion phase. Although accretion 
may decrease somewhat in Class I (or become episodic, 
see Sj8|), it appears unlikely that the average accretion 
rate drops by more than a factor of two from Class to 
Class I, based on the relative lifetimes of the two phases. 
Similar mean luminosities for the Class and Class I 
phases also argue against very high accretion rates in 
Class 0. 

Our derived Clas s lifetime is similar to the results of 
IVisser et al.l ()2002[ ) for a sample of Lynds dark clouds 
( teiass n ^ 2 X 1 0^ yr) , and to the recent findings 
of iHatchell et all ()2007ah that the Class hfetime in 
Perseus is similar to the Class I lifetime (2.5 — 6.7 x 10^ yr 
with 95% confidence). 

Our large, unbiased sample provides a distinct advan- 
tage over many other previous studies, which have nec- 
essarily relied on small samples or accumulated sources 
from a number of different surveys, wavel engths, and 
detect ion methods. For example, the Viss er et alJ 
(120021) Class samp le consists of 7 sources, and the 
lAndre fc Montmerld (|T994) lifetime for Ophiuch us is 
based on one Class object. IFroebrich et al.l ()2006D note 
that their 50 Class O/I sources are selected from a variety 



of surveys including NIR imaging of outflows, IRAS data, 
submillimeter and millimeter mapping, and radio contin- 
uum surveys, causing their source sample to be subject 
to strong selection effects. Within each cloud our surveys 
are very uniform, providing protostellar samples that are 
envelope mass limited and not biased by selection effects. 
For this reason, our estimated Class lifetimes should be 
more robust than most prev i ous me asurements. The re- 
cent work bv IHatchell et alJ ()2007al ) comparing SCUBA 
850 fim maps and Spitzer c2d data of Perseus, with 34 
Class sources, is a notable exception. 

Classifying sources based on T^'^;, calculated from pho- 
tometry corrected for extinction (see § I3.1.ip results in 
a somewhat smaller number of Class sources, and a 
slightly shorter Class lifetime. Table |7| gives the class 
statistics and corresponding Class lifetime derived us- 
ing the extinction corrected photometry. The three-cloud 
average Class lifetime is shorter than when using ob- 
served photometry, by approximately 35% (i^iass o = 
1.1 X 10^ yr). 

9.1. Cloud to Cloud Differences: Ophiuchus 

The results in Ophiuchus are strikingly different from 
Perseus and Serpens, with nearly 10 times more Class I 
than Class sources. Our results for Ophiuchus are sim- 
ilar to previously o bserved ratios of -/Vciass o/ -/Vciass i ^ 
1/10 in that cloud (I Andre fc Montmer lei [1991 . The de- 
rived Class lifetime, tciass o ~ 7 x lO'' yr, is still sub- 
stantially longer than the very shor t tciass o ~ lO'' yr 
found bv , Andre fc Montmerld (|1994t). however. 

There are two obvious, but conflicting, possible expla- 
nations for the smaller iVciass oZ-^Ciass i ratio in Ophi- 
uchus. First, the star fo rmation rate may b e temporally 
variable in Ophiuchus. Viss er et al.l ()2002l ) suggest that 
a burst of star formation approximately 10^ yr ago is re- 
sponsible for the large ratio of Class I to Class sources. 
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Fig. 15. — Evolution of the bolometric temperature with time, based on estimated hfetimes for four T^oj intervals from T^ai < 50 K to 
650 K. Tjoi ranges, inferred lifetimes, and uncertainties are taken from TableO and plotted as thick gray lines and symbols. The dashed line 
connects to the Class II lifetime at t = 2 X 10® yr and Tf,oi = 2800 K. We fit a simple function to Tt^i versus time: Tf,oi = 25 K+C(t/10^ yr)"; 
the best fit is for Thni oc t^-^. 



Alternatively, the Class phase may be much shorter in 
Ophiuchus than in the other two clouds due to higher ac- 
cretion rates at early times. Were this the case, however, 
we would expect the mean luminosity to be significantly 
higher for the Class sources than for the Class I sources, 
which is not observed. While a burst of star formation 
seems the more likely explanation, accretion rates, and 
thus lifetimes, could conceivably depend on cloud envi- 
ronmental factors such as mean density or turbulent ve- 
locity. 

On the other hand, the observed number of Class 
sources in Ophiuchus may be biased by the dearth of 
160 /im flux measurements in that cloud. As discussed 
in fAl bolometric temperatures are likely overestimated 
for sources without a 160 ^m fluxes, which may cause 
up to 50% of Class sources to be classified as Class I. 
The lack of 160 fim fluxes is especially problematic in 
Ophiuchus, where the majority of sources are either sat- 
urated at 160 /im or in regions of poor coverage, and may 
be partially responsible for the small number of observed 
Class sources. 

9.2. Limitations on the Class I Sample 

Recall that we only include sources in our Class I sam- 
ple if they are detected at 1.1 mm. This criteria excludes 
24 sources from the original candidate protostar sample 
in Perseus with 70 < Tbd < 650 K, 18 sources in Ser- 
pens, and 36 in Ophiuchus. The average spectrum of 



these "non-envelope Class I" sources is shown in Fig- 
ure [13] (Perseus and Serpens only) , and appears to be 
intermediate between the Class I and Class II averages. 
Furthermore, the mean Thoi values of the non-envelope 
Class Is (570 K in Perseus, 520 K in Serpens, and 500 K 
in Ophiuchus) confirm that they fall at the warm end of 
the Class I distribution. While it is possible that we are 
simply not detecting the envelopes of these sources (e.g., 
if Menv ^0.1 Mq or they are very diffuse), we suggest 
that the transition from Stage I (protostars that retain 
an envelope) to Stage II (pre-main sequence stars with 
no envelope) occurs closer to T^oi = 400 — 500 K than 
the standard T^oi = 650 K, for these data. 

When we add the "non-envelope Class I" sources to our 
Class I samples, dropping the requirement that they be 
detected at 1.1 mm, the Class to Class I ratio becomes 
0.5 in Perseus, 0.3 in Serpens, and 0.07 in Ophiuchus. 
This represents a 25-30% decrease i n the Class hfetime 
for each cloud. On the other hand, ICrapsi et al.l (j2008D 
suggest based on protostellar models that up to half of 
sources observationally classified as Class I may actually 
be Stage II source with disk inclinations > 65°. Our 
"non -envelope Cla s s I" o bjects may be just such sources; 
if the lCrapsi et al.l ()2008l ) picture is correct, the observed 
number of Class I sources likely represents an upper limit 
to the true number of Stage I protostars. 

9.3. Evolution of Tboi with Time 
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In Table [8] we give the number of sources in each cloud 
in our narrow Tboi bins ( "early Class 0" , "late Class 0" , 
etc.; ^jH]). The ratio of the total number in each bin, 
summed over the three clouds, to the total number of 
embedded protostars is used to derive a lifetime for each 
Tfioi range. The mean Tboi within that range is also given. 

With the exception of "early Class 0" , there is no ev- 
idence for a difference in the number of sources in each 
Tijoi bin, and little difference in the derived lifetimes. No- 
tably, 75% of the embedded phase lifetime has elapsed 
by the time the bolometric temperature reaches 300 K, 
and the mean T^oi in "Late Class I" (400 K) is skewed 
toward the lower temperature end of the bin. These 
features support our earlier suggestion that the divid- 
ing line between Stage I and Stage II is probably closer 
to Tboi = 400 - 500 K than to 650 K. Tboi = 50 K is 
reached at quite early times, which is not unexpected if 
the temperature scale starts at approximately 10 K with 
starless cores. 

With lifetime measurements in several Tboi intervals, 
together with our previous conclusion that the average 
accretion rate is approximately constant through the em- 
bedded phase, we can empirically characterize the evo- 
lution of Thoi with time. Figure [T5] plots the evolution 
of the bolometric temperature with time, based on val- 
ues from Table m The relationship is clearly non-linear; 
given the small number of observed points, we fit a simple 
function to T^oi as a function of time: 



T, 



bol 



25K + C 



t 



105 yr 



(7) 



where C is constant and we choose 25 K as the tempera- 
ture aX t = because that is the minimum observed T^oi 
in Class 0. The best fit is for C = 51 K and Ttoi cx 
indicating a fairly steep evolution with time. The depen- 
dence on time must flatten out significantly after the cm- 
bedded phase, however, to match the lifetime for Class II 
(dashed line in Figure fT5|). 

10. CONCLUSIONS 

Utilizing large-scale 1.1 mm surveys (jEnoch et al.l 
I2006H Young et al.ll2006HEnoch et afl2007f) together with 
Spitzer I RAC and MIPS m aps from the c2d Legacy 
program ()Evans et all l2003f l. we have constructed an 
unbiased census of deeply embedded protostars in the 
Perseus, Serpens, and Ophiuchus molecular clouds. Our 
sample includes a total of 39 Class sources and 89 
Class I sources, with approximate detection limits of 
Menv ^ 0.1 M© and L^oi > O.O5L0 for the envelope 
mass and bolometric luminosity, respectively. We also 
detect a few Class II and Herbig Ae/Be candidates at 
1.1 mm, most likely evidence of fairly massive proto- 
planetary disks. 

Bolometric luminosities, temperatures, and envelope 
masses are calculated for the candidate Class and 
Class I sources in each cloud. We compare protostel- 
lar classification methods, concluding that, for deeply 
embedded sources, the bolometric temperature T^oi is a 
better measure of evolutionary state than the near- to 
mid-infrared spectral index, am. 

We also explore classifying sources into "early Class 
0", "late Class I", etc., based on dividing them into nar- 
rower Thoi bins. Average observed spectr a in these bins 
are compared to model predictions from IWhitnev et al.l 



(|2003f) for "early Stage 0" , "late Stage 0" , etc. In a broad 
sense the Stage and Stage II models match the observed 
spectra fairly well. The agreement with both Stage I 
models is quite poor, however, as the models severely 
over-predict the near-IR and under-predict the mid-IR 
flux at this stage, displaying a double-peaked SED that 
is rarely observed. 

Observed source properties are compared to pro- 
tostellar evolutionary models using the bolometric 
temperature- luminosity (Lboi — Ti,oi) diagr am, the pro- 
tostellar equivalent of the H-R diagram (|Mvers et al.l 
Il998[ ) . Neither models wit h a constant mass accretion 
rate (lYoung fc Evansll20?)5l). nor those w ith an exponen- 
tially decreasing rate (jMvers et al .11 19981 ) fit the observed 
sources. 

In particular, there is a large population of low lu- 
minosity Class I sources that aggravate the previously 
note d "luminosity prob lem" for embedded protostars 
(e.g. lKenvon et al1ll990f ). We interpret this result as ev- 
idence for episodic accretion beginning at least by the 
Class I phase, and possibly earlier. More than 50% of 
Class I sources are inferred to have "sub-Shu" mass ac- 



cretion rates (M < 10 



M 



yr" 



corresponding to 



Lboi ^ lL0)i and approximately 20% have M < 10~ 
M© yr^^. To build up of order a solar mass in 5.4x 10^ yr, 
such sources must also have periods of "super-Shu" ac- 
cretion {M > 10~^ Lgyr^i). Few very high luminosity 
sources are observed (5%), suggesting that such rapid ac- 
cretion periods must be short lived. An important caveat 
to this analysis is that we may sometimes underestimate 
the luminosity due to missing 70 or 160 /zm photometry. 

Finally, the relative number of Class and Class I 
sources are used to estimate the lifetime of the Class 
phases. There are approximately half as many Class 
as Class I sources in the three cloud sample, implying an 
average Class lifetime of 1.7 ± 0.3 x 10^ yr (1.1 x 10^ 
yr when approximate extinction corrections are applied). 
This lifetime rules out drastic changes in the mass accre- 
tion rate from Class to Class I, particularly extremely 
rapid accretion in Class 0. In Ophiuchus the fraction of 
Class sources is much smaller. While this difference 
could be due to the lack of 160 /im flux measurements 
in Ophiuchus, it may be that either the Class phase is 
shorter in that cloud (0.7 x 10^ yr), or that a burst of 
star formation is responsible for the larg e population o f 
Class I objects (e.g. as suggested bv iVisser et al.l 120021 ). 

Altogether, the large variation in Class I luminosities, 
similar mean luminosities for the Class and Class I 
phases, and not dramatically different numbers of Class 
and Class I sources (at least in Perseus and Serpens) 
suggests a consistent picture of nearly constant average 
mass accretion rate through the entire embedded phase 
(Stage and Stage I), with highly variable episodic ac- 
cretion turning on by at least early Stage I, and possibly 
sooner. Understanding embedded protostellar structure 
and evolution well enough to reproduce the observations 
with detailed models of spectra and evolution presents 
an ongoing challenge. 
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TABLE 1 

Dkftnition of Clashes and Stacks 



Class/Stage 


definition 


Class 


Tioi < 70 K 


Class I 


70 K< Tboi < 650 K; 1.1 mm detection (Me„„ > O.IMq) 


Class II 


650 K< Ttoi < 2800 K 


Early Class 


n^i < 50 K 


Late Class 


50 K < Ttoi < 100 K 


Early Class I 


100 K < Tfco, < 300 K 


Late Class I 


300 K < Ttoi < 650 K 


non-envelope Class I 


70 K< Ttoi < 650 K; no 1.1 mm detection {Menv < O.IMq) 


Stage 


M-t < Menv 


Stage I 


iW, > Menv] Menv > O.lAf© 


Stage II 


circumstellar disk; M^nv < 0.1 Mpj 



TABLE 2 

BOLOMETRIC TEMPERATURES, LUMINOSITIES, AND ENVELOPE MASSES OF 
EMBEDDED PROTOSTARS IN PeRSEUS 



ID 




c2d name/position 


Tboi 


Lbol 


dm 


Mcnv 


Bolocam ID 


Other Names 






(SSTc2dJ...) 




(K) 


(Le) 






Mo) 








Class 


Per-emb 


1 


J034356.53-I-320052 


9 


24 (5) 


1.5 (0.4) 


0.5 (0.18) 


3.02 (0.08) 


Bolo 


103 


HH 211 


Per-emb 


2 


J033217.95-I-304947 


6 


25 (2) 


1.3 (0.2) 


1.07 (0.11) 


2.88 (0.08) 


Bolo 


66 


IRAS 03292-1-3039 


Per-emb 


3 


J032900.52-I-311200 


7 


30 (4) 


0.69 (0.15) 


2.16 


(0.1) 


0.29 (0.03) 


Bolo 


41 




Per-emb 


4 


J032839. 10-1-310601 


8 


31 (3) 


0.22 (0.03) 


1.68 


(0.06) 


0.36 (0.02) 


Bolo 


30 




Per-emb 


5 


J033120.96-I-304530 


2 


33 (1) 


1.2 (0.2) 


0.98 


(0.06) 


0.95 (0.04) 


Bolo 


65 


IRAS 032824-3035 


Per-emb 


6 


J033314.40-1-310710 


9 


36 (5) 


0.64 (0.01) 


2.22 


(0.06) 


1.04 (0.1) 


Bolo 


79^ 




Per-emb 


7 


J033032.68-I-302626 


5 


37 (4) 


0.15 (0.06) 


2.08 


(0.07) 


0.99 (0.08) 


Bolo 


62 




Per-emb 


8 


J034443.94+320136 


1 


40 (6) 


3.0 (0.5) 


0.96 


(0.07) 


0.63 (0.06) 


Bolo 


116=^ 




Per-emb 


9 


J032951. 82-1-313906 


1 


41 (2) 


0.46 (0.01) 


3.44 


(0.11) 


0.64 (0.04) 


Bolo 


59 


IRAS 03267-1-3128 


Per-emb 


10 


J033316.45-I-310652 


5 


47 (14) 


0.27 (0.66) 


1.73 


(0.09) 


1.57 (0.16) 


Bolo 


79=^ 




Per-emb 


11 


J034356.85+320304 


6 


49 (15) 


0.45 (0.85) 


1.37 


(0.09) 


1.78 (0.04) 


Bolo 


104 


IC 348-MMS 


Per-emb 


12 


J032910.50+311331 





51 (17) 


4.2 (4.4) 


2.58 


(0.08) 


7.75 (0.78) 


Bolo 


48a 


NGC 1333-IRAS 4A 


Per-emb 


13 


J032912.04+311301 


5 


54 (18) 


1.6 (2.0) 


0.55 


(0.06) 


3.66 (0.37) 


Bolo 


48^ 


NGC 1333-IRAS 4B 


Per-emb 


14 


J032913.52+311358 





54 (14) 


0.49 (0.32) 


2.41 


(0.13) 


0.5 (0.05) 


Bolo 


48^ 


NGC 1333-IRAS 4C 


Per-emb 


15 


J032904.05+311446 


6 


56 (16) 


0.53 (0.67) 


1.43 


(0.06) 


1.16 (0.12) 


Bolo 


46 




Per-emb 


16 


J034350.95+320324 


8 


56 (17) 


0.38 (0.62) 


1.51 


(0.06) 


1.15 (0.11) 


Bolo 


102=' 




Per-emb 


17 


J032739.09+301303 





59 (11) 


4.2 (0.1) 


2.68 


(0.06) 


0.51 (0.03) 


Bolo 


22 


RNO 15-FIR 


Per-emb 


18 


J032911. 25+311831 


3 


59 (12) 


2.8 (1.7) 


1.94 


(0.06) 


2.47 (0.07) 


Bolo 


49'' 




Per-emb 


19 


J032923.49+313329 


5 


60 (3) 


0.36 (0.05) 


1.51 


(0.06) 


0.64 (0.05) 


Bolo 


57 




Per-emb 


20 


J032743.23+301228 


8 


60 (14) 


1.7 (0.01) 


2.39 


(0.06) 


0.5 (0.03) 


Bolo 


23 


L1455-IRS 4 


Per-emb 


21 


J032910.68-I-311820 


5 


63 (12) 


2.8 (1.9) 


1.95 


(0.06) 


2.47 (0.07) 


Bolo 


49'' 




Per-emb 


22 


J032522.33-I-304514 





63 (11) 


1.7 (1.1) 


2.34 


(0.07) 


1.41 (0.14) 


Bolo 


5 


L 1448-IRS2 


Per-emb 


23 


J032917. 16-1-312746 


4 


66 (18) 


0.29 (0.33) 


1.75 


(0.06) 


0.5 (0.05) 


Bolo 


52 




Per-emb 


24 


J032845.30-I-310542 





67 (10) 


0.43 (0.01) 


1.11 


(0.05) 


0.19 (0.02) 


Bolo 


33* 
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TABLE 2 - Continued 



ID c2d name/position T^oi ^bol 

(SSTc2dJ...) (K) (Lq) 



Mcnv 

Mo) 



Bolocam ID 



Other Names 



Per-emb 25 J032637.46+301528.0 68 (12) 0.95 (0.02) 1.09 (0.05) 
Per-emb 26 J032538.82+304406.3 69 (7) 4.4 (1.5) 2.16 (0.06) 
Per-emb 27^^ J032855.56+311436.6 69 (1) 19.0 (0.4) 3.03 (0.12) 



0.5 (0.08) 
0.19) 
0.28) 



1.87 
2.81 



Class 1 


Per 


emb 


28 


J 034351 


02+320307 


9 


72 (25) 


0.28 (0.37) 


-0.28 (0.06) 





62 


(0.06) 


Bolo 


102" 




Per 


emb 


29 


J033317 


85+310932 





76 (20) 


1.8 (2.5) 


3.33 (0.06) 


3 


84 


(0.09) 


Bolo 


80 


Bl-c 


Per 


emb 


30 


J033327 


28+310710 


2 


78 (6) 


1.7 (0.01) 


1.93 (0.05) 





59 


(0.06) 


Bolo 


84 




Per 


emb 


31 


J032832 


55+311105 


2 


80 (13) 


0.16 (0.01) 


0.78 (0.05) 





19 


(0.02) 


Bolo 


25 




Per 


emb 


32 


J034402 


40+320204 


9 


84 (29) 


0.16 (0.31) 


1.53 (0.06) 





63 


(0.02) 


Bolo 


106 




Per 


emb 


33 


J032536 


48+304522 


3 


90 (18) 


4.3 (3.7) 


2.62 (0.07) 


4 


47 


(0.45) 


Bolo 


8 


L 1448-N 


Per 


emb 


34 


J033015 


12+302349 


2 


99 (13) 


1.6 (0.1) 


1.7 (0.05) 





46 


(0.07) 


Bolo 


60 


IRAS 03271+3013 


Per 


emb 


35 


J032837 


09+311330 


7 


103 (26) 


9.1 (0.3) 


2.35 (0.08) 





44 


(0.04) 


Bolo 


29 


NGC 1333-IRAS 1 


Per 


emb 


36 


J032857 


36+311415 


7 


106 (12) 


5.3 (1.0) 


1.6 (0.06) 





74 


(0.07) 


Bolo 


38=^ 




Per 


emb 


37 


J032918 


27+312320 





106 (21) 


0.07 (0.35) 


1.26 (0.06) 





44 


(0.04) 


Bolo 


54=^ 




Per 


emb 


38 


J033229 


18+310240 


9 


115 (21) 


0.54 (0.01) 


0.4 (0.05) 





45 


(0.08) 


Bolo 


68 




Per 


emb 


39 


J033313 


78+312005 


2 


125 (47) 


0.04 (0.08) 


1.44 (0.06) 





41 


(0.04) 


Bolo 


78 




Per 


emb 


40 


J033316 


66+310755 


2 


132 (25) 


1.5 (1.0) 


1.57 (0.05) 


1 


12 


(0.11) 






Bl-a 


Per 


emb 


41 


J033320 


30+310721 


3 


157 (72) 


0.17 (0.36) 


0.88 (0.05) 


1 


.5 (0.15) 


Bolo 


81 


Bl-b 


Per 


emb 


42 


J032539 


10+304358 





163 (51) 


0.68 (0.85) 


2.36 (0.06) 


1 


82 


(0.18) 


Bolo 


lO'' 




Per 


emb 


43 


J034202 


16+314802 


1 


176 (42) 


0.07 (0.06) 


1.47 (0.05) 





07 


(0.02) 








Per 


emb 


44c 


J032903 


76+311603 


7 


188 (9) 


32.5 (7.1) 


1.21 (0.08) 


3 


25 


(0.32) 


Bolo 


43 


SVS 13A 


Per 


emb 


45 


J033309 


57+310531 


2 


197 (93) 


0.05 (0.06) 


1.13 (0.05) 





21 


(0.02) 








Per 


emb 


46 


J032800 


40+300801 


3 


221 (7) 


0.3 (0.07) 


0.95 (0.05) 





17 


(0.02) 








Per 


emb 


47 


J032834 


50+310051 


1 


230 (17) 


1.2 (0.1) 


0.83 (0.05) 





12 


(0.02) 






IRAS 03254+3050 


Per 


emb 


48 


J032738 


23+301358 


8 


238 (14) 


0.87 (0.04) 


-0.19 (0.05) 





34 


(0.04) 


Bolo 


21 


L1455 FIR2 


Per 


emb 


49 


J032912 


94+311814 


4 


239 (68) 


1.1 (0.7) 


1.05 (0.05) 





77 


(0.08) 


Bolo 


49a 


ASR 30 


Per 


emb 


50^ 


J032907 


76+312157 


2 


254 (23) 


10.0 (3.0) 


2.18 (0.1) 


1 


62 


(0.16) 


Bolo 


47 




Per 


emb 


51 


J032834 


53+310705 


5 


263 (115) 


0.07 (0.1) 


0.54 (0.05) 





24 


(0.04) 


Bolo 


28 




Per 


emb 


52 


J032839 


72+311731 


9 


278 (119) 


0.16 (0.21) 


0.57 (0.05) 





49 


(0.05) 


Bolo 


31=" 




Per 


emb 


53 


J034741 


56+325143 


9 


287 (8) 


4.7 (0.9) 


0.78 (0.05) 


1 


67 


(0.08) 


Bolo 


122 


B5-IRS1 


Per 


emb 


54<i 


J032901 


57+312020 


7 


304 (63) 


6.2 (2.6) 


2.09 (0.08) 


1 


87 


(0.19) 


Bolo 


42 


NGC 1333-IRAS 6 


Per 


emb 


55 


J034443 


33+320131 


4 


309 (64) 


1.8 (0.8) 


0.5 (0.05) 





32 


(0.03) 


Bolo 


116=^ 


IRAS 03415+3152 


Per 


emb 


56 


J034705 


42+324308 


4 


312 (1) 


0.54 (0.09) 


0.36 (0.05) 





14 


(0.02) 






IRAS 03439+3233 


Per 


emb 


57 


J032903 


33+312314 


6 


313 (200) 


0.09 (0.45) 


1.13 (0.05) 





11 


(0.02) 








Per 


emb 


58 


J032858 


44+312217 


4 


322 (88) 


0.63 (0.47) 


0.83 (0.05) 





71 


(0.07) 


Bolo 


40 




Per 


emb 


59 


J032835 


04+302009 


9 


341 (179) 


0.04 (0.06) 


0.15 (0.05) 





14 


(0.02) 


Bolo 


27 




Per 


emb 


60 


J032920 


07+312407 


5 


363 (240) 


0.28 (1.05) 


0.42 (0.05) 





52 


(0.05) 


Bolo 


54=* 




Per 


emb 


61 


J034421 


33+315932 


6 


371 (107) 


0.24 (0.16) 


0.21 (0.05) 





17 


(0.02) 


Bolo 


113 




Per 


emb 


62 


J034412 


98+320135 


4 


378 (29) 


1.8 (0.4) 


0.37 (0.05) 





12 


(0.02) 








Per 


emb 


63 


J032843 


28+311733 





436 (9) 


1.9 (0.4) 


0.36 (0.05) 





25 


(0.02) 


Bolo 


31=^ 




Per 


emb 


64 


J033312 


85+312124 


1 


438 (8) 


3.2 (0.6) 


0.41 (0.05) 





53 


(0.07) 


Bolo 


76 




Per 


emb 


65 


J032856 


31+312227 


8 


440 (191) 


0.16 (0.16) 


-0.14 (0.05) 





25 


(0.02) 








Per 


emb 


66 


J034345 


15+320358 


6 


542 (110) 


0.69 (0.22) 


-0.22 (0.05) 





24 


(0.02) 









Bolo 18 
Bolo lO'' 
Bolo 38=^ 



L 1448-C 

NGC 1333-IRAS 2A 



^ NOTE — The envelope mass Menv is calculated from the 1.1 mm flux assuming a dust temperature To = 15 K. Bolocam IDs are 
from Papers I— III, and indicate that the protostellar source is within 1.0 core FWHM of a given Bolocam core. If no ID is listed, the 
source was not originally identified as a distinct 1.1 mm core, but there is detectable (> 3a-) 1.1 mm emission at the position of the source. 
Uncertianties for Ti,o;, Li,o;, and Menv are given in parentheses; uncertainties for Tjoj and I/i,oi are taken to be the difference of values 
calculated by the midpoint and prismoidal methods (see CA)| . Uncertainties in Menv are from the photometric errors in the total flux only. 
Absolute uncertainties in the mass may be a factor of 2 or more, from uncertainties in d, K\rnm and T15, but relative values should be much 
more accurate. 

" More than one protostar is associated with this Bolocam core or source is in a crowded region; mass is based on flux measured in a 
30 — 40" diameter aperature. 

^ More than one protostar is associated with this Bolocam core; mass is based on dividing Bolocam core flux equally (source separation 
is less than 30"). 

Saturated 24/.tm flux has been replaced with IRAS 25/im flux for the calculation of Z/joj and T(,oj . 

Saturated at either 24^m or 70^m so Ltoj is likely an underestimate, but the IRAS beam is too confused to measure reliable 25/im or 
60/im fluxes. 
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TABLE 3 

BOLOMETRIC TEMPERATURES, LUMINOSITIES, AND ENVELOPE MASSES OF COLD PROTOSTAHS IN SERPENS 



ID c2d name/position T^oi ^bol ctlR Mcnv Bolocam ID Other Names 

(SSTc2dJ...) (K) (Lq) (Mq) 



Class 



Ser-emb 1 


J182909.24-I-003132.3 


39 (2) 


1.6 (0.1) 


2.27 (0.08) 


1.16 


(0.02) 


Bolo 15 




Ser-emb 2 


J182952.44-I-003611.7 


42 (10) 


1.0 (0.1) 


0.77 (0.05) 


0.54 


(0.03) 


Bolo 24=^ 




Ser-emb 3 


J182854.84-I-002952.5 


51 (12) 


2.6 (0.1) 


1.91 (0.06) 


1.54 


(0.15) 


Bolo 8^ 




Ser-emb 4 


J183000.72-I-011301.4 


54 (16) 


1.2 (1.6) 


1.68 (0.08) 


2.56 


(0.03) 


Bolo 28=" 




Ser-emb 5 


J182854.84-I-001832.6 


56 (12) 


0.16 (0.1) 


0.9 (0.06) 


0.24 


(0.02) 


Bolo 7=* 




Ser-emb 6 


J182949.56-I-011521.9 


56 (12) 


11.0 (6.0) 


2.65 (0.07) 


7.98 


(0.07) 


Bolo 23^ 


FIRSl 


Ser-emb 7 


J182854.12-I-002929.9 


58 (13) 


3.1 (0.1) 


1.36 (0.06) 


1.67 


(0.17) 


Bolo 8^ 




Ser-emb 8'' 


J182948.12-f011644.9 


58 (16) 


2.1 (2.4) 


1.37 (0.06) 


3.72 


(0.37) 


Bolo 22=^ 


S68N 


Ser-emb 9 


J182855.92-I-002944.7 


66 (21) 


0.34 (0.63) 


1.89 (0.06) 


1.19 


(0.12) 


Bolo 8=^ 





Class I 



Ser- 










5 




1 HQ Id Cl^) 


1 '^'^ (f) OR"! 


0.29 


(0.03) 


Bolo 


3a 


IRAS 1862-1-0050 


Ser- 


-emb 


11 


J182906, 


.72-f003034, 


,3 


77 (12) 


1.9 (1.2) 


1.66 (0.06) 


1.39 


(0.14) 


Bolo 


14a 






Ser- 


-emb 


12*' 


J182952, 


,08-1-011547, 


,8 


85 (9) 


2.6 (1.2) 


1.54 (0.06) 


1.23 


(0.12) 


Bolo 


23=^ 


SMM 10 IR 




Ser- 


-emb 


13 


J182902, 


,04-1-003120, 


,6 


86 (31) 


0.05 (0.09) 


0.24 (0.06) 


0.17 


(0.02) 










Ser- 


-emb 


14 


J183005, 


,40-f004104. 


,5 


100 (32) 


0.08 (0.1) 


1.3 (0.07) 


0.19 


(0.02) 










Ser- 


-emb 


15 


J182954, 


,24-1-003601, 


,3 


101 (43) 


0.17 (0.27) 


-0.18 (0.06) 


0.48 


(0.05) 


Bolo 


24a 






Ser- 


-emb 


16'' 


J182844, 


,76-1-005125, 


,7 


110 (26) 


0.04 (0.6) 


1.07 (0.05) 


0.24 


(0.02) 


Bolo 


3=^ 






Ser- 


-emb 


17 


J182906, 


,36-1-003043, 


,2 


117 (21) 


1.5 (1.3) 


1.7 (0.06) 


1.38 


(0.14) 


Bolo 


14a 






Ser- 


-emb 


18 


J182952, 


,80-f011456. 


,0 


120 (15) 


1.2 (0.8) 


1.49 (0.05) 


0.74 


(0.07) 


Bolo 


23=^ 






Ser- 


-emb 


19 


J183000, 


,00-f011311. 


,6 


129 (29) 


2.0 (1.9) 


2.57 (0.06) 


2.56 


(0.03) 


Bolo 


28=^ 






Ser- 


-emb 


20*^ 


J182949, 


,20-f011619. 


,8 


131 (3) 


4.4 (2.1) 


3.8 (0.11) 


1.53 


(0.15) 


Bolo 


22'' 






Ser- 


-emb 


21 


J182951, 


,00-1-011640, 


,6 


141 (37) 


1.7 (1.5) 


1.04 (0.05) 


1.84 


(0.18) 


Bolo 


22'' 






Ser- 


-emb 


22 


J182957, 


,48-1-011300, 


,5 


157 (39) 


2.6 (2.2) 


1.01 (0.05) 


2.35 


(0.24) 


Bolo 


25" 






Ser- 


-emb 


23 


J182957, 


,84-1-011251, 


,4 


190 (23) 


3.1 (1.5) 


0.67 (0.05) 


1.06 


(0.11) 


Bolo 


25" 






Ser- 


-emb 


24 


J183000, 


,00-1-011158, 


,9 


249 (64) 


1.8 (1.2) 


1.15 (0.05) 


1.14 


(0.11) 


Bolo 


29 






Ser- 


-emb 


25 


J182851, 


,24-1-001927, 


,3 


250 (106) 


0.06 (0.07) 


0.45 (0.05) 


0.16 


(0.02) 


Bolo 


7" 






Ser- 


-emb 


26 


J182958, 


,92-1-011426, 


,2 


303 (139) 


0.62 (0.73) 


0.46 (0.05) 


2.0 


(0.2) 


Bolo 


26" 






Ser- 


-emb 


27'' 


J182956, 


,76-1-011446, 


,5 


322 (25) 


5.8 (1.7) 


0.3 (0.05) 


0.69 


(0.07) 










Ser- 


-emb 


28'' 


J182844, 


,04-f005337. 


,9 


360 (244) 


0.19 (1.5) 


0.45 (0.05) 


0.26 


(0.03) 


Bolo 


2 






Ser- 


-emb 


29 


J183002, 


,88-1-011228, 


,2 


374 (101) 


2.2 (1.3) 


0.21 (0.05) 


1.32 


(0.13) 


Bolo 


28" 






Ser- 


-emb 


30 


J182957, 


,84-1-011405, 


,7 


560 (85) 


13.8 (4.2) 


0.28 (0.09) 


1.9 (0.19) 


Bolo 


26" 


IRAS 18274H 


-0112 


Ser- 


-emb 


31 


J182931, 


,92-1-011842, 


,9 


430 (15) 


5.3 (1.0) 


0.26 (0.05) 


2.12 (0.07) 


Bolo 


20 






Ser- 


-emb 


32 


J183005, 


,76-f003931. 


,7 


431 (207) 


0.05 (0.07) 


-0.39 (0.05) 


0.2 (0.02) 


Bolo 


32 






Ser- 


-emb 


33<= 


J182916, 


,08-1-001822, 


,7 


437 (10) 


1.7 (0.2) 


-0.07 (0.05) 


0.54 (0.04) 


Bolo 


17 


IRAS 18267H 


-0016 


Ser- 


-emb 


34 


J182902. 


,76-1-003009, 


,5 


631 (410) 


0.12 (0.28) 


-0.14 (0.05) 


0.22 (0.02) 


Bolo 


13 







Note. — Columns are as in Table[2] and uncertainties for Tfjoj , Lf^^l: ^-nd Menv are given in parentheses. 
^ More than one protostar is associated with this Bolocam core or source is in a crowded region; mass is based on flux measured in a 
30 — 40" diameter aperature.^ Saturated at either 24fim or 70/j-m so Lfj^i is likely an underestimate, but the IRAS beam is too confused 
to measure reliable 25/im or 60/im fluxes.'^ Missing 70/j.m flux has been replaced with IRAS 60//m flux for the calculation of L^oi and 
'^boi-'^ Source lies outside the 70/im map.*^ Saturated 24/im flux has been replaced with IRAS 25/xm flux for the calculation of Lh^i and 

Tbal. 
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TABLE 4 

BOLOMETRIC TEMPERATURES, LUMINOSITIES, AND ENVELOPE MASSES OF COLD PROTOSTARS IN OPHIUCHUS 



ID c2d name/position T^bol ^hol ctlR Mcnv Bolocam ID Other Names 

(SSTc2dJ...) (K) (Lq) (Mq) 



Class 



Oph-cmb 1 


J162821. 72-243623. 4 


35 (2) 


0.25 (0.01) 


1.23 (0.06) 


0.51 


(0.04) 


Bolo 26 


IRAS 16253-2429 


Oph-emb 2=* 


J163222.56-242831.8 


54 (8) 


6.9 (1.2) 


5.03 (0.2) 


1.16 


(0.02) 


Bolo 36 


IRAS 16293-2422B 


Oph-emb 3 


J162625.80-242428.8 


57 (10) 


0.41 (0.15) 


1.65 (0.07) 


0.24 


(0.02) 


Bolo 9 


VLA 1623 (?) 








01 


ass 1 














Oph-emb 


4 


J163136.84-240419.9 


77 (25) 


0.18 (0.05) 


-0.27 (0.05) 


0.03 


(0.01) 








Oph-emb 


5 


J162721. 96-242727.7 


87 (28) 


0.08 (0.08) 


-0.05 (0.05) 


0.12 


(0.01) 


Bolo 


17 




V.7 Uli eiliU 


g 


J162705. 40-243629. 5 


106 (14) 


0.15 (0.09) 


1.27 (0.05) 


0.09 


(0.01) 


Bolo 


14b 


LFAM26/GY197 


Oph-emb 


7 


J163151. 96-245725. 9 


124 (39) 


0.1 (0.13) 


0.82 (0.05) 


0.05 


(0.01) 


Bolo 


33 






gc,d 


11 62621 48-242304 2 


1 33 f34") 


18.3 (1.4) 


1 46 Co 09 ~1 


0.22 


(0.02) 


Bolo 


5 


OSS 30 


Oph-emb 


9 




1 (aqS 


0.1 (0.3) 


U . / \\J . \JO J 


0.53 


(0.05) 


Bolo 


8 


GY 30 


Oph-emb 


10 


1163200 96-245642 7 




2.7 (0.3) 


1 (0 07) 


0.05 


(0.01) 


Bolo 


35 


L1689S/TRS67 


Oph-emb 


11 


J162717.64-242856.2 


190 (1) 


0.52 (0.06) 


0.25 (0.05) 


0.07 


(0.01) 








Oph-emb 


12 


1162724 48-244103 1 


191 (30 ~l 


33 (0 16~1 


1 01 (n 05") 


0.12 


(0.01) 


Bolo 


1?^ 




Oph-emb 


13<= 


J162728. 08-243933. 4 


193 (31) 


7.2 (0.5) 


2.29 (0.11) 


0.11 


(0.01) 






IRS 44/GY269 ) 


Oph-emb 


14c 


J162727.00-244050.5 


211 (7) 


3.8 (0.7) 


1.17 (0.08) 


0.36 


(0.04) 


Bolo 


18*= 


IRS 43/GY265 


Oph-emb 


15 


J163152.32-245536.1 


240 (42) 


0.13 (0.07) 


1.07 (0.05) 


0.08 


(0.01) 








Oph-emb 




J162709.36-243718.4 


257 (60) 


17.9 (0.6) 


1.69 (0.15) 


0.07 


(0.01) 


Bolo 


14b 


EL 29/GY214 


Oph-emb 


17 


J163135.76-240129.2 


290 (18) 


1.5 (0.1) 


0.14 (0.05) 


0.14 


(0.02) 


Bolo 


30 


IRAS 16295-2355 


Oph-emb 


18 


J162857.72-244054.8 


304 (96) 


0.03 (0.03) 


0.67 (0.05) 


0.04 


(0.01) 








Oph-emb 


19 


J162728.44-242720.8 


309 (85) 


0.47 (0.28) 


-0.03 (0.05) 


0.29 


(0.03) 


Bolo 


20*= 


IRS45/GY273 


Oph-emb 


20 


J162706.84-243815.0 


310 (23) 


0.6 (0.16) 


0.61 (0.05) 


0.05 


(0.01) 


Bolo 


13 


WL17/GY205 


Oph-emb 


21<: 


J162702.16-243727.1 


331 (2) 


5.2 (1.3) 


1.53 (0.11) 


0.08 


(0.01) 






WL16/GY182 


Oph-emb 


22 


J162640.56-242714.4 


372 (160) 


0.06 (0.09) 


0.45 (0.05) 


0.06 


(0.01) 






GY91/CRBR42 


Oph-emb 


23 


J162648.48-242838.6 


428 (84) 


0.12 (0.05) 


0.02 (0.05) 


0.03 


(0.01) 






WL2/GY128 


Oph-emb 


24 


J162737.08-244237.8 


438 (87) 


0.13 (0.06) 


0.13 (0.05) 


0.04 


(0.01) 






GY301 


Oph-emb 


25 


J163143.68-245524.6 


496 (68) 


0.28 (0.09) 


0.23 (0.05) 


0.04 


(0.01) 








Oph-emb 


26 


J162730.24-242743.2 


528 (99) 


0.93 (0.27) 


-0.12 (0.05) 


0.19 


(0.02) 


Bolo 


20'' 


IRS47/GY279 


Oph-emb 


27 


J162739.96-244314.8 


551 (54) 


0.74 (0.15) 


-0.15 (0.05) 


0.04 


(0.01) 






IRS51/GY315 


Oph-emb 


28 


J162721. 60-244143.0 


595 (57) 


1.2 (0.2) 


-0.03 (0.05) 


0.04 


(0.01) 






IRS42/GY252 



Note. — Columns arc as in Tablc[2] and uncertainties for T^o/ , L^oi, and Menv are given in parentheses. 
^ Oph-cmb 2 is saturated at 160/J-m; when the lOO/J-m IRAS flux is included the resulting L^oj is 16±6Lq, consistent with the luminosity 
calculated by Wa lker et al.| jl986D when scaled to a distance of 125 pc.^ More than one protostar is associated with this Bolocam core 
or source is in a crowded region; mass is based on flux measured in a 30 — 40" diameter aperature.*^ Saturated 24/im flux has been 
replaced with IRAS 25/J-m flux for the calculation of Lboi and Thoi-'^ Saturated TOfim flux has been replaced with IRAS 60/j,m flux for the 
calculation of L^o/ and Tj^^i ■ 



TABLE 5 

Numbers of Protostars by Class and the Class lifetime 



Cloud 


NciassO 


NciassI 


Nciasso/NciassI 


TClassO 












(yr) 


Perseus 


27 (27) 


39 


(41) 


0.70 (0.66) 


2.2 ±0.4 X 10^ 


Serpens 


9 (13) 


25 


(23) 


0.4 (0.6) 




Ophiuchus 


3 (4) 


25 


(24) 


0.12 (0.17) 


0.7l:o 5 X 10^ 


Three cloud sample 


39 (44) 


89 


(88) 


0.44 (0.50) 


1.72 ±0.25 X 10^ 



Note. — Numbers of Class and Class I sources are based on T^oi classifications fTablcfTl 
including the requirement that Class I sources be detected at 1.1 mm. Numbers in parentheses 
indicate how the results change if we utilize a different method for calculating Ttoi (prismoidal 
versus midpoint integration, see ^A\ . The Class lifetime, TclassO^ assumes th at the entire 
embe dded phase lasts for 5.4 x 10'^ years {Temb — '''ciassO-\-'''ciassi — 5.4 x 10^ vr: |Evans et all 
I2OOI1 . and Nciasso/iNc lassQ + ^Ciassi) — "^classO I Temb- The lifetime is calculated using the 
average of the two Nciasso/ ^Classi ratios given, and the uncertainty is from the difference 
between them added in quadrature with the y/N uncertainties from counting statistic. 
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TABLE 6 

Class II sources detected at 1.1 mm in Perseus, Serpens, and Ophiuchus 



c2d name/position Tboi Ljjoi 
(SSTc2dJ...) (K) (Lq) 



am 



(Mo) 



Bolocam ID'' 



other names 



J032747.66+301204.6 662 (15) TEO 

J034400.00+320154.1 700 (17) 1.60 

J032900.00+312146.8 856 (26) 0.33 

J034109.13+314438.0 873 (30) 1.70 

J032856.64+311835.6 993 (22) 0.42 

J034548.27+322411.8 1096 (23) 5.00 

J032917.66+312245.1 1297 (22) 1.50 



-0.09 (0.05) 
-0.34 (0.05) 
-0.83 (0.05) 
-0.72 (0.05) 
-0.78 (0.05) 
-0.72 (0.05) 
-0.77 (0.05) 



0.06) 
0.04) 
0.09) 
0.02) 
0.02) 
0.02) 
0.02) 



RNO 15^= (1) 

LkHa 353 (2) 
IRAS 03380-1-3135 
ASR 120 (3) 
LkHa 330 (2) 
LkHa 270 (2) 



0.05) 
0.03) 
0.01) 
0.06) 
0.01) 
0.09) 
0.02) 



0.82 
0.40 
0.90 
0.16 
0.10 
0.10 
0.19 



Per-Bolo 24 
Per-Bolo 106* 
Per-Bolo 40* 
Per-Bolo 90 
Per-Bolo 39* 
Per-Bolo 120 
Per-Bolo 54* 



J182900.00+003003.1 806 (5) OJW 

J183006.12-I-004233.9 871 (2) 3.60 

J182900.96-f002931.6 892 (2) 1.20 

J182901.32-f002933.0 892 (6) 0.36 

J182901. 68-1-002954.6 1158 (3) 1.80 

J182901. 68-1-002946. 5 1212 (5) 0.70 



Serpens 



-0.61 (0.05) 
-0.42 (0.05) 
-0.42 (0.05) 
-0.54 (0.05) 
-1.04 (0.05) 
-0.65 (0.05) 



0.03) 
0.01) 
0.02) 
0.02) 
0.02) 
0.02) 



Ser-Bolo 13* 
Ser-Bolo 33 
Ser-BoIo 13* 
Ser-Bolo 13* 
Ser-Bolo 13* 
Ser-Bolo 13* 



0.01) 
0.01) 
0.01) 
0.01) 
0.01) 
0.01) 



0.35 
0.17 
0.23 
0.23 
0.23 
0.25 



Ophiuchus 



CoKu Ser G3 (4) 
CoKu Ser G4 (4) 



J163151. 96-245615.7 
J162658.56-244531.6 
J162624.00-241613.4 
J163945.36-240203.8 
J162623.28-242100.7 



690 (71) 
824 (14) 
958 (35) 
970 (30) 
1080 (100) 



0.3) 

0.1) 

0.2) 

0.09) 

0.4) 



-0.17 (0.05) 
-0.45 (0.05) 
-0.71 (0.05) 
-0.73 (0.05) 
-0.62 (0.05) 



0.01) 
0.01) 
0.01) 
0.01) 
0.01) 



1.2 
1.3 
1.9 
0.84 
2.0 



0.05 
0.03 
0.04 
0.03 
0.15 



GY 186 (5) 



Oph-Bolo 7* GY 20 (5) 



Note, — References: (2) IHerbig fc Kameswarall972l : Ol IAspin et al,ll994l : ('4 HCohen fc Kuhai979l : (Sl IGreene fc Yound 



ne 1.1 mm mass (Mimm) is calculated for an optically thin disk and a dust temperature of Tx) — 15 K. Our assumption 
for the dust opacity may be incorrect for typical Class II disks, however, where grain growth is a possibility fe.g., lN'atta et al.l 
120071) Bolocam IDs with "*"s indicate that the source is in a confused region at 1.1 mm, and the association of the Spitzer 
object with the given 1.1 mm emission is not secure.*' Typically classified as a Class I source, RNO 15 has a bolometric 
temperature just above the Class I cutoff in our sample. The spectral index of RNO 15 is cxjr — —0.09, placing it in the 
"fiat spectrum" category, or close to the Class I boundary. 



TABLE 7 

Effect of Extinction Corrections on Class Statistics 



Cloud 


NciassO 


NciassI 


NciassO /NciassI 


TClassO 










(yr) 


Perseus 


21 


42 


0.5 


1.5 X 10^ 


Serpens 


4 


26 


0.15 


0.6 X 10^ 


Ophiuchus 


3 


20 


0.15 


0.6 X 10^ 


Three cloud sample 


29 


88 


0.33 


1.1 X 10^ 



Note. — Same as Table \E\ but for classifications based on T^^i , the extinction- 
corrected bolometric temperature. In addition, an embedded phase lifet ime of 4.4 X 
10^ y r is assumed, derived based on extinction-corrected classifications dEvans et al.| 
f200l) . 



TABLE 8 

Lifetimes for narrow T^^i bins 



Phase 


Ttoi range 
(K) 


N(Per) 


N(Ser) 


N(Oph) 


N(total)/ 
N(cmb,tot) 


Mean Uoi 
(K) 


Lifetime 

(y) 


Early Class 


T < 50 


11 


2 


1 


14/128 = 0.11 


40 


0.6 X 10^ 


Late Class 


50 < T < 100 


23 


11 


4 


38/128 = 0.30 


70 


1.6 X 10^ 


Early Class I 


100 < T < 300 


19 


12 


12 


43/128 = 0.34 


180 


1.8 X 10-^ 


Late Class I 


300 < T < 650 


13 


9 


11 


33/128 = 0.26 


400 


1.4 X 10^ 



Note. — For each Tjoj bin, N(total) is the total number of sources in all three clouds; N(emb,tot) 
is the total number of embedded protostars in the three cloud sample (128). Only sources detected at 
1.1 mm are included. The mean Tjoj of protostars in a given Tjoj range is also given. Lifetimes are 
calculated from N(total)/N(emb, tot) X 5.4x 10^ yr. Cloud-to-cloud variations are typically 0.2 — 0.3 X 10^ 



Protostars in Serpens, Perseus, and Ophiuchus 



29 



APPENDIX 

CALCULATING THE BOLOMETRIC LUMINOSITY AND TEMPERATURE 

Determination of the bolometric temperature and luminosity of any given source can depend strongly on the method 
used to approximate the integrations over frequency in equations[TH31 As 5*^ is sampled at a finite number of frequencies 
(usually 6 to 10 here), the SED must be interpolated over the intermediate frequencies. Here we use two different 
methods; the first (midpoint) method utilizes a simple linear interpolation for the midpoint flux, while the second 
(prismoidal) method calculates a color temperature for each pair of flux points, and uses a modified blackbody based 
on that color temperature to estimate the midpoint flux density. In both cases the SED is extrapolated from the 
longest wavelength flux using Su , and flux upper limits are removed from the fit (i.e., we interpolate over them). 
Examples of the midpoint and prismoidal interpolations are shown in Figure [TB] for two sources with Class 0-like (left) 
and Class I-like (right) SEDs. 
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Fig. 16. — Examples of the interpolated midpoints used for the integration of the SED in the midpoint (solid line) and prismoidal (dashed 
line) methods. Observed data points are overlaid as diamonds. A typical Class source is shown on the left, and a Class I source on 
the right; both sources shown are lacking data at 160 /xm. After the midpoint is calculated, the integration for both methods is done in 
linear vSv space. The prismoidal method provides a good approximation of the dust emission peak for cold sources (left), but probably 
overestimates the long- wavelength flux for warmer, flatter SED sources (right), resulting in a bias toward lower Ti,o;. 



Figure [17] compares the results of using the midpoint and prismoidal methods to calculate Lboi and Tboi for the 
embedded protostar sample in Perseus. The midpoint and prismoidal T^oi values agree fairly well, usually to within 
20%, with no strong systematic bias. There is some tendency for the prismoidal T^oi to be lower for colder sources. 
On the other hand, L^oi calculated by the prismoidal method is generally higher than the midpoint L^oi- 

To further test this behavior, we plot in Figure [TH] the fractional errors due solely to finite sampling, for L^oi and 
Tboi calculated from input blackbody spectra. Blackbody sources with input T^oi from 5-5000 K are sampled at the 
observed wavelengths (A =1.25, 1.65, 2.17, 3.6, 4.5, 5.8, 8, 24, 70, 160, 1100 /xm), and the bolometric temperature and 
luminosity calculated by the midpoint (left) and prismoidal (right) methods. Both Lboi and Tboi have typical errors 
of approximately 20%, and the measured Tboi is consistently over-estimated by approximately 20%. The midpoint 
Tboi is more variable but somewhat more accurate than the prismoidal Tboi^ while Lboi is somewhat less accurate. Of 
course, the blackbody spectra tested here are not necessarily representative of the more complicated observed SEDs. 
Based on Figures [TTHTSI and on examinations of the midpoint and prismoidal fits to observed SEDs, from which we 
find that the prismoidal method often provides a poor fit to "flatter" SEDs at long wavelengths (Figure [TB]) , we use 
the midpoint method throughout this paper. Differences between the midpoint and prismoidal methods are often used 
as a measure of uncertainty. 

In general, measured 160 /im fluxes are uncertain by up to a factor of 2, due to unquantified saturation effects 
and calibration uncertainties, and in many cases we are unable to measure a 160 /xm flux at all. In addition, the 
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Midpoint Tt,„| (K) Midpont Ti,^, (K) 

Fig. 17. — Comparison of the midpoint and prismoidal methods for calculating Li,ol and T(,o;, for sources in Perseus. Left: Midpoint and 
prismoidal Ttoi (black asterisks) and Lt,ol (triangles) values are plotted for all candidate embedded protostars in Perseus. Right: Fractional 
difference between the midpoint and prismoidal methods, plotted as (midpoint-prismoidal)/prismoidal. The two methods generally agree 
to within 20% for Tjoj but Ljoj calculated by the prismoidal method is higher than the midpoint L^qj by as much as 100%. There is some 
tendency for the prismoidal Ti,o; to be lower for colder sources. 

"point-source" 160 /ini fluxes of extended sources will likely be underestimated, because they are measured with a 
PSF fit rather than aperture photometry. Missing or severely underestimated fluxes at 160 /xm can significantly affect 
both Lboi and Ti,oi for cold sources that peak near 100 /xm. As we are primarily interested in the coldest sources, we 
attempt here to quantify the effects of missing or underestimated 160 /xm flux. One way of approaching this problem 
is to examine sources for which we ha ve additional information near the peak of the SED, e.g., from IRAS or SHARC 
350 ^m observations ()Wu et al.ll2007f) . 

We first look at two example sources, one cold (IRAS 03282+3039; Tbo; ~ 33 K), and one warmer {Thoi ~ 163 K). 
IRAS 03282+3039 is isolated, so we can use the IRAS fluxes without worrying about confusion due to the large IRAS 
beam. We calculate the "true" Ltoi and Tfcoz by including the IRAS 60 and 100 /um fluxes, a SHARC II 350 /im 
flux, and a MIPS 160 /im flux measured in a large aperture, which is almost two times higher than the PSF-fit flux. 
Using the underestimated point source 160 /im flux, and no longer including IRAS or SHARC II data, causes an 
underestimate of Ltoi by 35% and an overestimate oiTboi by 15% compared to the "true" values. Omitting the 160 fiia 
point altogether, as would be appropriate for saturated sources, results in an underestimate of Lhoi by 7% and an 
overestimate of T^oi by 9%. 

Thus we conclude that while the 160 /im point is important for characterization of the SED of embedded objects, 
our integration method can interpolate over a missing 160 fim flux to estimate Lhoi and Thoi to within 20% for cold 
sources. Severely underestimated 160 /im flux densities will cause larger errors of up to 50%. Using a similar procedure 
for the warmer source (Tboi — 163 K), we find that omitting the 160 /im fiux results in errors in the opposite sense 
compared to the colder source: an overestimate of Lboi by 28% and an underestimate of T^oi by 18%. 

For a more general result, we compare Thoi calculated with the 160 /im fiux included in the SED to the value found 
by omitting the 160 /im point, for all sou rces in Perseus wit h a measured 160 /im flux (Figure [191 left panel). For this 
plot, we have adopted classiflcations from lChen et al.l ()1995[ ). as discussed in §31 Although Tftoj calculated without the 
160 /tm flux is always overestimated compared to Tboi calculated with the 160 /im flux, it is never a large enough effect 
to shift the source classification of objects with Tboi ^ 100 K (e.g., from Class I to Class II). This is not the case for 
Class sources, however, as can be seen in the right panel of Figure [191 where the scale has been adjusted to highlight 
the lowest Tboi values. Five sources that have Tboi < 70 K (Class 0) when including the 160 /im fiux are shifted to 
650 < Tboi > 70 K (Class I) when the 160 /im fiux is omitted. If a published 350 fim fiux is available (diamonds), the 
errors in Tboi resulting from excluding the 160 /im fiux are almost completely eliminated. 

Based on Figures [l7liT9l and the above examples, we estimate overall uncertainties for measured Lboi and Tboi values 
of 20 — 50%, depending on whether or not a 160 /tm fiux is available. If we are unable to measure a 160 /tm flux, Tboi 
will almost certainly be an overestimate for very cold sources, which may affect our classification of the most deeply 
embedded protostars. 
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Fig. 18. — Characterization of sampling errors for T[,o; and L^o; ealeulated using the midpoint (left) and prismoidal (right) methods. The 
fractional difference between the input (Tboj ,Li,o;) of blackbody spectra and the measured values when SEDs are sampled at the observed 
wavelengths is plotted as a function of the input Tij^i- Both Lj,oi T^qj have typical errors of approximately 20%, and the measured 
Ti,o! is consistently over-estimated by approximately 20%. Source SEDs, which are considerably more complex than the blackbodies tested 
here, may behave differently. 




T^Qi with 1 60/Am Tj,Q| with 1 60 /im 

Fig. 19. — Characterization of the systematic errors introduced into the measured T[,o; when a 160 /im flux is not available. Left: For all 
sources in Perseus with a reliable 160 /im flux, we calculate Tf,ol with and without the 160 /im point (asterisks). Tf,ol always overestimated 
when the 160 ^m flux is missing, but it does not change the classification of Class I sources for these data. Right: Similar, but with the 
scale adjusted to highlight lower Ti^oi sources. Not having a 160 /xm point does change the classification of approximately half of the Class 
{I'bol < 70 K) sources. Errors axe substantially reduced when a flux at 350 ixm. is available (diamonds). 



